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We present an analysis of the structure and properties of chromatic polynomials 
P(G p t,rn,q) of one-parameter and multi-parameter families of planar triangulation 
graphs G ptj rn, where fh = (mi, ...,m p ) is a vector of integer parameters. We use these 
to study the ratio of \P(G pt ^,T + 1)| to the Tutte upper bound (r — l) n ~ 5 , where 
t = (1 + >/5)/2 and n is the number of vertices in G ptt rn- In particular, we calculate 
limiting values of this ratio as n — > oo for various families of planar triangulations. 
We also use our calculations to study zeros of these chromatic polynomials. We 
study a large class of families G vt ,m with p = 1 and p = 2 and show that these have 
a structure of the form P(G p t^ m , q) = c G iAj 1 + c G?rf) 2A™ + c G?rf) 3A™ for p = 1, where 
Ax = q-2, A 2 = q-3, and A 3 = -1, and P(G pt ^q) = £? 1=1 £* =1 c Gpt ^ A™ 1 A™ 2 
for p = 2. We derive properties of the coefficients c ^ and show that P{G pt ^q) 
has a real chromatic zero that approaches (1/2) (3 + y/b ) as one or more of the 
mi — > oo. The generalization to p > 3 is given. Further, we present a one-parameter 
family of planar triangulations with real zeros that approach 3 from below as m — > oo. 
Implications for the ground-state entropy of the Potts antiferromagnet are discussed. 



I. INTRODUCTION 



In this paper we present exact results on the structure and properties of chromatic poly- 
nomials P{Gpt,mi q) °f one-parameter and multi-parameter families of planar triangulation 
(pt) graphs G ptt rn, where m = (m-i, m p ) is a vector of integer parameters. Our results 
substantially generalize our previous study in [lj. In standard notation we let G = (V,E) 
be a graph with vertex and edge sets V and E, and denote the number of vertices and edges 
as n = n(G) = \V\ and e(G) = \E\, respectively. The resultant set of faces is denoted 
F(G), with cardinality f(G) = \F(G)\. We recall the definitions of a planar graph G as one 
that can be drawn in a plane without any crossing edges, and a triangulation as a graph all 
of whose faces are triangles. For an arbitrary graph G, the chromatic polynomial P(G, q) 
enumerates the number of ways of associating q colors with the vertices of G, subject to 
the constraint that adjacent vertices have different colors (called a proper g-coloring of G) 
[2]. The minimum number of colors for a proper g-coloring of a graph G is the chromatic 
number of G, denoted x(G). Without loss of generality, we restrict here to graphs G that 
are connected, have no multiple edges, and have no loops (where a loop is defined as an edge 
that connects a vertex back to itself). Multiple edges are also excluded by our restriction 
here to triangulation graphs, since they produce faces that are not triangles. An important 
identity is P(G,q) = Z P ^ F (G,q,0), where Z P AF(G,q,T) denotes the partition function of 
the Potts antiferromagnet (PAF) on the graph G at temperature T. Because of this identity, 
properties of chromatic polynomials are of interest both for mathematical graph theory and 
for statistical physics. A particularly significant feature of the Potts antiferromagnet is the 
fact that it generically exhibits nonzero ground-state (i.e., zero-temperature) entropy per 
vertex for sufficiently large q on a given graph. 

The chromatic polynomial of a graph G may be computed via the deletion-contraction 
relation P(G, q) = P(G — e, q) — P(G/e, q), where G — e denotes G with an edge e deleted and 
G/e denotes the graph obtained by deleting e and identifying the vertices that it connected. 
From this follows the cluster formula P(G,q) = J2 g > C g(~ l) 6 *- '- 1 q k{ - G '\ where G' = (V,E') 
with E' C E and k(G') denotes the number of connected components of G'. Via this relation, 
one can generalize q from positive integers to real and complex numbers, as is necessary when 
analyzing zeros of P(G, q), called chromatic zeros. For a general graph G or family of graphs 
Gtft, the calculation of the chromatic polynomial takes an exponentially long time. Planar 
triangulations comprise a class that is particularly amenable to analysis, as we shall show. 
Since a triangulation graph G t (whether planar or not) contains at least one triangle, P(G t , q) 
always contains the factor P(K 3 , q) = q(q — l)(q — 2). (Here, K n is the complete graph with 
n vertices, defined as the graph such that each vertex is adjacent to every other vertex by 
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an edge.) 

An interesting upper bound on an evaluation of a chromatic polynomial of a planar 
triangulation was derived by Tutte [3] (see also [I], [3]), namely 

0< \P(G pt7 r+l)\ <U(n(G pt )) , (1.1) 

where r = (1 + Vo)/2 is the golden ratio and 

U(n) = r 5 - n = (r - l) n " 5 . (1.2) 

As in PQ, it is natural to define the ratio 

r(G ) ~ ' P ^ r + 1 )l fl 3) 

r{Gpt) = U(n(G pt )) ' (L3) 

which is bounded above by 1. In this paper we shall present calculations of this ratio, and its 
limit as n — > oo, for a number of different families of planar triangulations. In p] we showed 
that if p = 1, then, with m = m, if P(G p t, m , q) involves only a single power of a polynomial, 



(Ac p t) m as in (2.8) below, it follows that r(G p t :7n ) approaches zero exponentially rapidly as 



m — > oo. In [T] we also constructed one-parameter families G pt) m for which P(G p t, m} q) is 



a sum of powers of certain terms A» with i = 1,2,3, given below in (3.1), then r(G p t t m) 
may approach a finite nonzero constant as m — > oo. We generalize these results here to the 
families G pt ^ with p > 2. We also exhibit a p = 1 family, denoted F m , for which P(F m , q) is 
a sum of nonpolynomial A^jS, i = 1,2,3, and show that for this family, r(F m ) decreases to 
zero (exponentially rapidly) asm-) oo. It should be noted that the Tutte upper bound is 
satisfied as an equality by the triangle, K3, but for planar triangulations with higher n, the 
upper bound is realized strict inequality. 

Part of our work concerns zeros of chromatic polynomials P(G p t.m, q) for families of planar 
triangulations G pt> rn, i-e., chromatic zeros of these graphs. An interesting aspect of this study 
relates to an empirical observation made by Tutte in connection with his upper bound, 
namely that for a planar triangulation G pt , P(G pt , q) typically has a real zero close to q — 
r + 1 = 2.6180339887.... This observation has been somewhat mysterious over the years, for 



at least two reasons. First, Tutte's upper bound (1.1) does not imply that a P(G pt ,q) need 
have a zero near to r+1. Second, it is known that for an arbitrary (loopless) graph G, P(G, q) 
cannot vanish exactly atg = r + l = (3 + \/5 )/2. We recall the proof. If P(G, (3 + \/5 )/2) 
were zero, then, since P(G, q) is a polynomial in q, it would have the factor [q — (3 + \fh )/2] . 
But since P(G, q) has rational coefficients (actually, integer coefficients, but all we use here 
is the property that they are rational), it would therefore also have to contain a factor 
involving the algebraic conjugate root, namely [q — (3 — \/5 )/2]. However, this would imply 
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that P(G, q) would also vanish at q = (3 — \fh )/2 = 0.381966..., but this is impossible, since 
this point lies in an interval (0,1) where P(G,q) cannot vanish J6j[7]. As part of our analysis 
here, we shed some light on this mystery by constructing families of planar triangulation 
graphs Gpt^m with p = 1 and p = 2 for which the chromatic polynomials have the respective 
forms P(G p t 7 m,q) = c Gpt ,i^T'^ c a pt ,^2 iJrC a pt ,^3 > " for p — 1 (with mi = m), where Ai = q — 2, 
A 2 = q-3, and A 3 = -1, and P(G pttA , q) = E? 1= i Ei 2 =i ^.Wa^ 1 ^ for P = 2 - We derive 
properties of the coefficients c * and show that P(G p trn, q) has a real chromatic zero that 

approaches (1/2) (3 + \fh ) as one or more of the — > oo. We give the generalization of 
this result to p > 3. Further, we construct a one-parameter family of planar triangulations 
of this type with real zeros that approach 3 from below as m — > oo. 



II. GENERAL PROPERTIES OF ONE-PARAMETER FAMILIES OF PLANAR 

TRIANGULATIONS 



We have constructed and studied various one-parameter families of planar triangulations 
Gpt^m that can be built up in an interative (recursive) manner. (Here and below, for families 
where fh is one-dimensional, we set m\ = m to simplify the notation.) In this section we de- 
rive general properties of the chromatic polynomials of these families of planar triangulations. 
For our families, the number of vertices is linearly related to m, 



where a and (3 are constants that depend on the type of family. We recall the Euler relation 
|V(Cr)| — | | + |F(G)| = xe = 2 for a graph G embedded on a surface of genus 0, such as 
the plane, where xe is the Euler characteristic. In general, for a planar graph each of whose 
faces has p sides, n(G), e(G), and f{G) satisfy the relations e(G) = p(n(G) — 2)/(p — 2) 
and f{G) = 2(n(G) — 2)/{p — 2). For the case of interest here, namely planar triangulation 
graphs, where each face is a triangle, it follows that 



A. 



General 



n(G pt ^ m ) =am + (3 , 



(2.1) 




(2.2) 



and 




(2.3) 



so that e{G pt ) = (3/2) f(G pt ). 
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In our present study we will make use of several results that we derived in [T|. First, since 
U(n) — > as n — > oo and since m is proportional to n, it follows that for these families of 
planar triangulations, 

lim P(G pt , m ,r + l) = . (2.4) 



Second, given the upper bound (1.1) and the fact that U(n) approaches zero exponentially 
fast as n — > oo, it follows that 

P(G p t t m, t + 1) approaches zero exponentially fast as m — > oo. (2.5) 

We recall two definitions that apply to any graph: (i) the degree d(vi) of a vertex Vi G V 
is the number of edges that connect to it, and (ii) a fc-regular graph is a graph for which 
all vertices have degree k. Since a triangulation graph Gt is not, in general, /c-regular, it 
is useful to define an effective vertex degree in the limit \V\ — > oo. For this purpose, we 
introduce, as in our earlier work, the notation {G} for the formal limit n — > oo of a family 
of graphs G. We define 



d eff ({G}) = lim 



2\E\ 



|V|-wo \V\ 



(2.6) 



\V\ ' 

where for a given G, denotes the number of vertices with degree di and n{G) = \V\. 



Substituting (2.2) in (2.6), we obtain 

de ff ({G pt }) = 6 . (2.7) 



We will use this below, in Sect. XIX 



B. Families with P{G p t, m ,<l) Consisting of a Power of a Single Polynomial 

There are several ways of constructing one-parameter families of planar triangulations. 
One method that we have used is the following, which produces families for which the 
chromatic polynomial involves a single power of a polynomial in q. Start with a basic graph 
G p t } i, drawn in the usual explicitly planar manner. The outer edges of this graph clearly 
form a triangle, K 3 . Next pick an interior triangle in G ptt i and place a copy of G pty \ in this 
triangle so that the intersection of the resultant graph with the original G pt ^ is the triangle 
chosen. Denote this as G pt ^- Continuing in this manner, one constructs G pt , m with m > 3. 
The chromatic polynomial P(G ptj2 , q) is calculated from P(G ptt i, q) by using the s = 3 special 
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case of the complete-graph intersection theorem. This theorem states that if for two graphs 
G and H (which are not necessarily planar or triangulations), the intersection G(lH = K s for 
some a, then P(G UH,q) = P(G, q)P(H, q)/P(K a , q). (Note that P(K„ q) = UjZl(l ~ j)-) 
It follows that for planar triangulations formed in this interative manner, the chromatic 
polynomial has the form 

P(G pt , m ,q) = c Gpt (X Gpt r . (2.8) 

where the coefficient cc pt and the term Xc pt are polynomials in q that do not depend on m. 
Here and below, it is implicitly understood that m > m m i n , where m m i n is the minimal value 
of m for which the family G p t jm is well defined. Since G p t contains at least one triangle, K%, 
the coefficient c Gpt contains (and may be equal to) P(K 3 , q) = q(q — — 2). The chromatic 
number of G pt may be 3 or 4. In the case of a planar triangulation which is a strip of the 
triangular lattice of length m vertices with cylindrical boundary conditions, to be discussed 
below, an alternate and equivalent way to construct the (m + l)'th member of the family is 
simply to add a layer of vertices to the strip at one end. 



C. Families with P(G p t, m , q) Consisting of Powers of Several Functions 

We have also devised methods to obtain families of planar triangulations G pt ^ with the 
property that the chromatic polynomial is a sum of more than one power of a function of q. 
We begin with the simplest case, p — 1, i.e., one-parameter families and then discuss families 
with p > 2. For one-parameter families, we find the general structure 

Jmax 

P(G pt>m ,q) = J2 c o pt A X ^) m , (2-9) 

where m > m m i n and the c G t j and \G pt ,j are certain coefficients and functions depending 
on q but not on m. Here we use the label G pt to refer to the general family of planar 
triangulations G pt ^ m . We will describe these methods below. Parenthetically, we recall that 



the form (2.9) is a general one for one-parameter recursive families of graphs, whether or 



not they are planar triangulations [TT], [8]. For (2.9) evaluated at a given value q = qo, 
as m — > oo, and hence n — > oo, the behavior of P(G p t, m ,q) is controlled by which \G pt ,j 
is dominant at q — qo, i.e., which of these has the largest magnitude \^G pt ,j(qo)\- F° r our 
purposes, a go of major interest is r + 1, since the Tutte upper bound (1.1) applies for this 
value. We denote the \G pt ,j that is dominant at q = t + 1 as \G pt ,dom- Clearly, if P(G pt)m , q) 
involves only a single power, as in (2.8), then \ Gvt ,dom = ^G pt - 



As in earlier works [8T[TU]. it can be convenient to obtain the chromatic polynomials 
P(G p t,m, q) via a Taylor series expansion, in an auxiliary variable x, of a generating function 
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T(G p t, q, x). Below, we will have occasion to use this method for the family F m (see (18.4)). 



Both the form (2.9) and the expression via a generating function are equivalent to the 



property that P(G pttm ,q) satisfies a recursion relation, for m > j max + m r 



Jmax 



P(G pt:Tn , q) + 2J b Gpuj P(G pt , m -ji q) — 

3=1 



(2.10) 



where the b Gvtij , j = 1, ...j max , are given by 



Jmax 



Jmax 



3=1 3=1 



x 



Thus, 



Jmax 

ba pt ,i = ~ 2j ^G pt ,j , 
3=1 

Jmax 

i=i, fc=i, j^fe 



and so forth, up to 



pt Umax 



Jmax 

-iy max n 

i=l 



(2.11) 

(2.12) 
(2.13) 

(2.14) 



D. Asymptotic Behavior as m — > 00 

In pQ we discussed the asymptotic behavior of the chromatic polynomials as m — > 00. In 



both the cases of Eq. (2.8) and (2.9), a single power [XG pt j] m dominates the sum as m — > 00 



For a member of a one-parameter family of planar triangulations, G p t, m , we use the notation 

(2.15) 



r(G p t : m) for the ratio (1.3), and we define 

r(G pt!00 ) = lim r(G ptm ) . 
We define the (real, non-negative) constant a Gpt as p] 

,1/n _ |-V?p f ,rf m(T + l)| 1/a 



a Gpt = lim [r(G pt>m )\ 



t - 1 



(2.16) 



We showed in pQ that if j max = 1, then a Gpt < 1 and hence for the classes of G ptt m under 
consideration, (i) r(G p4i00 ) = and (ii) r{G pt> m) decreases toward zero exponentially rapidly 
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as a function of m and n as m — > oo. Note that this does not imply that P(G p t, m , q) has a 
zero that approaches g = r + 1 as m, n — >■ oo. 

For one-parameter families of planar triangulation graphs G pt>m where P(G p4jm , q) has the 



form (2.9) with j max > 2, a consequence of the Tutte upper bound (j X . 1 p is that as m — > oo, 
any contribution c g ^j(Xg v 



equal in magnitude to (r — l) n 



in (2.9), when evaluated at q = r + 1, must be less than or 



'herefore, for a given j in this case, either the coefficient 



c G j vanishes for q = t + 1 or, if this coefficient does not vanish at q — r + 1, then, taking 



into account the relation (2.1), it follows that 



I A 



•G P t,j 



1/a 



T- 1 



< 1 



at q = t + 1 V j . 



(2.17) 



If this inequality is realized as an equality, then r(G ptj00 ) is a nonzero constant, which nec- 
essarily lies in the interval (0, 1), so that ac pt = 1. For the families G pt}m for which m and 



n are linearly related, as specified in (2.1), this type of behavior occurs if and only if, when 
P{G p t,miQ) is evaluated at q — r + 1, the (necessarily) dominant Xq ptj j (with nonvanishing 
coefficient ca pt j), is equal to r — 1 in magnitude, i.e., \\G pt 4om\ = t — latg = r + l. This is 
true, in particular, if this Xc pt ,j = q — 2. In the structural form (3.16) below, we shall label 
this as the j = 1 term. 

It is a general property that if a graph G contains a complete graph K p as a subgraph, 
then P(G,q) contains the factor P(K p ,q). In particular, a triangulation graph, whether 
planar or not, has the factor P(K 3 , q) and a planar triangulation may also contain a K±. 
(However, by Kuratowski's Theorem, a planar graph may not contain a K p with p > 5.) 
Thus, for a planar triangulation G pt , P(G pt , q) = for q = 0, 1, 2. If G pi D K±, then 
P(G p t,q) also vanishes at g = 3. If P(G pt ,q) has the form of a single power, given as Eq. 
(2.8 ), then these factors are explicit. If, however, P(G pt , q) has the form of a sum of j max > 2 
powers [X Gpt ,j] m , then the conditions that P(G pt ,q) vanish at q — 0, 1, 2 imply relations 
between the various terms. Moreover, the condition that P(G ptjm , t + 1) obeys the Tutte 



upper bound (1.1 ) also implies conditions on the structure of this chromatic polynomial. We 



derive these next. 
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III. PROPERTIES OF A CLASS OF G pt , m WITH P(G pt , m ,q) HAVING j m ax = 3 

AND CERTAIN X Gpt>j 

A. Structure of Coefficients c G j in -P(G p t jm ,g) 

For a large class of one-parameter families of planar triangulations that we have con- 



structed and studied, for which P(G pt ,m, q) has the form (2.9), we find that (i) j max = 3 and 
(ii) the A G j = A, with j = 1, 2, 3 have the form 

pt 

Xx = q-2, A 2 = g-3, A 3 = -1 . (3.1) 

For this class of planar triangulations, we can derive some general results concerning the 
functional form of the coefficients c G j (where we will often suppress the subscript pt on G pt 



where the meaning is obvious). Using the general form (2.9) with j max = 3 and these X/s 



we can derive the following identities. The fact that P(G pt , 0) = implies that 

cg,i (-2) m + c G , 2 (-3) m + c G , 3 (-l) m = • (3.2) 

where for ease of notation we suppress the subscript pt on G pt here and in related equations. 
Since this equation must hold for arbitrary m (understood implicitly to be an integer in 
the range m > m min , where m min is the minimal value for which the family G pt , m is well 
defined), it implies that cqj = for all j. Hence, for these families, 

CQj contains the factor q for j — 1, 2, 3 . (3.3) 

The evaluation P(G p t, 1) = reads 

c Gil (-l) m + c G , 2 (-2) m + c Gi3 (-l) m = at 5 = 1. (3.4) 

Since this equation must hold for arbitrary m > m m i n , it implies two conditions on the 
evaluation of the coefficients at q — 1, namely 

c g,2 = and c G) i + c G3 = at 9 = 1. (3.5) 



In particular, (3.5) implies that 

c Gj2 contains the factor q — 1 . (3.6) 
The evaluation P(G pt , 2) = reads 

c G ,iO m + [c G , 2 + c G , 3 ](-l) m = at 5 = 2. (3.7) 
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Since this equation holds for arbitrary m > m m i n} it implies that 

cg,2 + c G ,3 = at q = 2 . 



(3i 



0, (3.7) and (3.8) together would also imply that cq,i = 



If a family G p t. m which has m m i 
at q = 2. 

Continuing with P(G p t }m , q) of the form ( 2.9 ) with (3.1), we next analyze the evaluation of 
P(G pt , m ,q) at g = r + 1, viz., P(G pt , m , r + 1). Since r-1 = 0.61803... and r-2 = -0.381966 
are smaller than unity in magnitude, the first two terms in P(G pt , m , r + 1) vanish exponen- 
tially rapidly as m increases. As regards the ratio r(G pttm ), as m increases, the contribution 
of the first term to this upper bound approaches a constant, while the contribution of the 
second term vanishes exponentially rapidly. Given the relation (2.1), the Tutte upper bound 
also vanishes exponentially rapidly as a function of m. Therefore, in order for P(G p t im , r + 1) 

(3.9) 

3 + y/E 



to satisfy the Tutte upper bound (1.1), it is necessary and sufficient that 

c G ,3 = at q = r + 1 . 

This means that 



cg,3 contains the factor q — (- 



(3.10) 



Given that a chromatic polynomial has rational (actually integer) coefficients as a polynomial 
in q, this means that cg,3 must also contain a factor involving the algebraically conjugate 
root, i.e., 

'3 - V5 \ 



Cg j3 contains the factor q — 
Combining these, we derive the result that 

Cq 3 contains the factor q 2 



3g + l 



(3.11) 



(3.12) 



Having proved these results, it is thus convenient to extract the factors explicitly and 
define 

(3.13) 



«G,2 = 



cga 



Q 

CG,2 



and 



«G,3 



q(q-l) 

CG,3 



(3.14) 



(3.15) 



q(q 2 — 3g + l) 

For this class of planar triangulation graphs G p t, m , we thus have the general structural 
formula 

P(G pt , m ,q) = qU G , 1 (q-2) m + K G , 2 (q-l)(q-3) m 
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+ K G , 3 (g 2 -3g+l)(-i; 



(3.16) 



where m and n are related by (2.1). We observe that the form (3.16) satisfies the general 
results that we derived above for the evaluation at q = r + 1. Thus, if P(G p t, m , q) has this 
form (3.16) with (2.1) and a — 1, then 



r(G. 



pi,oo y 



[9 K G,1 



g=r+l 



and hence 



a(G. 



The conditions on the coefficients cgj's evaluated at q = 1 and g 



(3.17) 
(3.18) 

2 that we have derived, 



(3.5), together with the definitions (3.13 )-(3.15 ), are equivalent to the following relations 



and 



«g,i = k g,3 at q = 1 



«g,2 = k Gi3 at g = 2 



(3.19) 



(3.20) 



For certain families of planar triangulations G p t im = G m , the chromatic number x(G m ) 
is 3 for even m and 4 for odd m or vice versa. In these cases, we can also derive another 
relation between the coefficients. Thus, if x(G m ) = 3 for even m and x(G m ) = 4 for odd m, 
then kg,i = «g,3 at g = 3. On the other hand, if x(G m ) = 3 for odd m and x(G m ) = 4 for 
even m, then re^i = — kg,3 at g = 3. In the case of families G m for which x(G m ) = 4 for all 
m, we have 

kg,i + ^g,3 (-l) m = at g = 3 if = 4 , (3.21) 



which implies 



kg, i = «G,3 = at g = 3 if x{G m ) = 4 



(3.22) 



B. Properties of Real Chromatic Zeros 

Here we derive some properties of chromatic zeros of planar triangulation graphs G 



pt,m 

for which the chromatic polynomial has the form (2.9). It is appropriate first to review some 
relevant properties of chromatic zeros of general graphs and planar triangulation graphs. For 
a general graph G, it is elementary that there are no negative chromatic zeros and that there 
are no chromatic zeros in the intervals (0,1) P, [7j. The property that (0,1) is a zero-free 
interval for the chromatic polynomial implies that q = r + 1 cannot be a chromatic zero for 
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any graph G, as noted above (independent of whether it is a planar triangulation or not). 
Another interval that has been proved to be free of chromatic zeros is (1, 32/27] [T21 [T3] . For 
an arbitrary graph G, let us denote the total number of subgraphs H C G that are triangles 
as N t , and let n = n(G) and e = e(G). It has been proved [J4] that for an arbitrary graph 
G with n > 4 vertices, if N t < u(G), where 

then P(G, q) has complex zeros. 

Specializing now to planar triangulation graphs, it has been proved that G p t has no 
chromatic zeros in the interval (2, q w ) [15], [16], where q w is the unique real zero of 

\ TC = q 3 - 9q 2 + 29q - 32 , (3.24) 

i.e., 

q w = 3 - l 12 ( 9 + )] 1/3 + 4 [ 12 ( 9 + v ^77 )]-i/3 = 2 .546602.. (3.25) 

We remark that q w occurs as a chromatic zero of some planar triangulation graphs, in par- 
ticular, the family comprised of cylindrical sections of the triangulation lattice with L y = 3, 
or equivalently iterated octahedra. In 1992 Woodall conjectured that a planar triangula- 
tion has no chromatic zeros in the interval (g m ,3), where q m = 2.6778146.. is the unique 
real zero of g 3 — 9q 2 + 30g — 35, [IB] , but later he gave counterexamples to his conjecture 
involving one-parameter families of planar triangulations each of which has a real zero that 
approaches 3 from below as this parameter goes to infinity [17] . We also note that for a 



planar triangulation, substituting (2.2) into (3.23), 



u(G) = (3n 2 ^ n 2) 5) for G = Gpt . (3.26) 

Here we present some further results on chromatic zeros of planar triangulations. First, 
if P(G ptt m,q) has the form (2.8) involving only a single power of a Ac pt , then its zeros are 



fixed, independent of m, and hence although it typically has a zero close to r + 1, this zero 
does not move as a function of m. However, if P(G p t, m ,q) has the multi-term form (3.16) 



with j max = 3 and the Aj's in (3.1), then it necessarily has a zero in the interval [q W} 3) that 
approaches r + 1 as m — > oo. The proof of this is as follows. Let us assume that q is a 
real number in this interval [^,3). In the limit as m — > oo, the first two terms, which are 
proportional to (q — 2) m and (q — 3) m , respectively, vanish (exponentially fast), so that 

P(G pt , m , q)~e m + q(q 2 - 3q + l)K Gpu3 (-l) m , (3.27) 
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where e m denotes the contribution of these first two terms. If /?gl ( ,3 happens to vanish at 
q = t + 1, then the result follows, since e m — > as m — > oo. If ^a vt ;i 7^ at g = r + 1, then 
consider the limit as q — > r + 1, where we can write 



■l) m P(G pt ,m,q) _ (-1) ,„ + g^_3g + i 



,2 



(r + l)K GpuZ \ q=T+1 (r + l)K Gpu3 \ q=T+1 

= 5 m + q 2 - 3q + 1 . (3.28) 
To show that P(G pt , q) has a zero that approaches g = r + 1 as m - > 00, we use the fact 



that P(G pt! m, l) is a continuous function of q and solve (3.28) for q, subject to the condition 
that q G [g w ,3), obtaining a consistent result with 

1 



q= 2 



3 + a/5 - 4<L 



(3.29) 



which approaches g = r+lasm— y 00. Note that the other zero at q = (l/2)(3 — \/5 — 4S m ), 
is irrelevant because we assumed at the outset that q is in the interval [q w , 3) and this other 
zero is outside this interval; in the vicinity of this other zero, the analysis does not apply 
because the terms proportional to (q — 2) m and (q — 3) m do not vanish as m — > 00. 



In 1.1 we have exhibited two one-parameter families of planar triangulation graphs with 
this property, namely B m and H m . We construct and analyze several more families of this 
type here. Thus, for these families, we have provided an understanding of why P(G pt! m, <?) 
has a chromatic zero near to r + 1 and, furthermore, have proved that this zero approaches 
r+ 1 as m — > 00. From the derivation above, it is evident that our result requires, for a given 
family G p t, m , that m be sufficiently large. As is illustrated from numerical results presented 
below, for specific families of planar triangulations that we have studied, P(G p t, m , q) has a 
real zero reasonably close to r + 1 even for moderate values of m. 

Our second result follows immediately from this analysis. With the same assumptions, we 
have observed that in the limit m — > 00, P(G pttm , <?) has a real zero in the interval q G [q w , 3) 
if and only if c G ti 3 has a real zero in this interval, q G [q w , 3). We know that there is at least 
one such zero, namely the one arising from the factor q 2 — 3q + 1 in c G Therefore, in the 
limit m — > 00, P(G p t, m , q) has another real zero in the interval q G [q w , 3) in addition to the 
one approaching r + 1 if and only if kg,3 has a real zero in this interval q G [q w , 3). We will 
present several applications of these results below. 

We remark that, with the same assumptions as above, 

lim \P(G pttm ,q)\ = c Gpt , 3 • (3.30) 

m— 5>oo 

Note that the limit linim^oo P(G p4 m , q) itself does not exist, because the term A™ = (— l) m 
factor has no limit as m — > 00. 
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In P, we investigated the question of whether for a planar triangulation graph G pt it is 
true that the chromatic zero of G pt nearest to r + 1 is always real. We exhibited an example, 
with a graph we denoted Gqm,i, which is, to our knowledge, the first case for which the zero 
closest to r + 1 is not real but instead the zeros closest to r + 1 form a complex-conjugate 
pair. Our result is in agreement with a previous observation by Woodall that this graph has 
no real zero near to r + 1 [18]. Since P(GcM ! m,Q) has the form of (2.8) with a single A, as 
m increases, its zeros are fixed and just increase in multiplicity, in contrast to the motion of 
chromatic zeros for families G p t >m whose chromatic polynomials are of the form (2.9) with 

jmax ^* 1 • 

The value of q where the Tutte upper bound applies, namely g = r + l = (3 + \fh )/2 is 
also a member of a sequence of numbers related to roots of unity, namely the Tutte-Beraha 
numbers, q r . Thus, for a root of unity of the form z r = e m ^ r , one defines q r = (z r + z*) 2 = 
4cos 2 (7r/r). One has r + 1 = q$. Parenthetically, we note that chromatic zeros have been 
studied for sections of triangular lattices with various boundary conditions that are not 
planar triangulations, either because they have at least one face that is not a triangle or 
because they are not planar (e.g., have toroidal or Klein-bottle boundary conditions). We 
refer the reader to [1] for references to some of these papers; here, in view of our focus on 
the Tutte upper bound (1.1) we restrict to planar triangulations. 



C. Properties of Complex Chromatic Zeros 



As before, we consider a one-parameter of planar triangulation graphs G p t t m such that 
P{G p t,m,q) has the form (3.16). Here we give a general determination of the continuous 



accumulation set B of chromatic zeros of P{G p t- m -, q) in the complex q plane in the limit 
m — > oo (and, hence, owing to (2.1), n — > oo). If the zeros form a discrete set (some with 



multiplicities that go to infinity as n — > oo), then this locus is null. As in earlier work, 
we denote the formal limit of the family G m as m — > oo as {G}. In general, the locus B 
may or may not intersect the real q axis. If it does, the maximal point where it intersects 



the real axis is denoted q c ({G}). For a P(G,q) of the form (2.9), the curves comprising 
the locus B are determined as the solutions of the equality in magnitude of the dominant 
Xoptj, in accordance with general results for recursive functions [IH]. These curves extend 
infinitely far from the origin if and only if such an equality can be satisfied as \q\ — > oo, as 
was discussed in [5] .[2D] -[22]. 

We now consider families of planar triangulations G ptt m whose chromatic polynomials 
have the form (3.16). For these families, first, because Ai = q — 2 is dominant for large 
q and is equal in magnitude to A 3 = — 1 at q = 3, it follows that q c ({G P t}) = 3. Second, 
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as a consequence of the fact that the equality \q — 2| = \q — 3| holds for the infinite line 
Re(g) = 5/2, part of the boundary B extends infinitely far away from the origin in the q 
plane, i.e., passes through the origin of the 1/q plane. As is evident for specific families G pt > m , 
as m and hence n go to infinity, the degree of one or more vertices also goes to infinity, so 
the fact that the magnitudes of zeros diverge is in accord with the upper bound \q\ < bA(G) 
obtained in [23J (with b ~ 7.96) and strengthened slightly in [23] (with b ~ 6.91), where 
A(G) is the maximal degree of any vertex in G. Note, however, that the property that the 
degree of a vertex diverges as n — > oo does not, by itself, imply that B passes through the 
origin of the 1/q plane. This is clear from the n — > oo limit of wheel graphs, for which the 
central vertex v cent . has degree d{v cent ) — > oo, but the locus B has bounded support in the q 
plane [HI 125] - Continuing with our analysis of families of planar triangulations G pty m whose 
chromatic polynomials have the form (3.16), the B separates the q plane into three regions, 
which we denote as Rj, j = 1, 2, 3. These are defined as follows: 

5 

i?i : Re(q) > - and \q - 2| > 1 , (3.31) 



2 
5 

R 2 : Re(q) < - and \q - 3| > 1 , (3.32) 



and 

i? 3 :|g-2|<l and \q - 3| < 1 . (3.33) 
The boundaries between these regions are thus the two circular arcs 

B(R\, R 3 ) : q = 2 + e ie , -| < 9 < | (3.34) 

and 

In 4tt 

B(R 2} R 3 ) : q = 3 + , — < <f> < — , (3.35) 
together with the semi-infinite vertical line segments 

B(R 1 ,R 2 ) = {q}: Re(g) = ^ and |Im(g)| > ^ . (3.36) 
These meet at the triple points 

Qt, Qt = g • ^ ' 

A specific example of a locus B that extends infinitely far from the origin in the q plane was 
studied in [HI [26] for the family B m . Here we have generalized the result to all families of 
planar triangulations whose chromatic polynomials have the form (3.16). 

We have constructed and studied a two-parameter family of planar triangulations D mi rn2 . 
By keeping one of the two indices mi or m 2 fixed, we have obtained a number of one- 
parameter families and have analyzed the chromatic polynomials for these. We have also 
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considered the special cases where one allows both mi and rri2 to vary, such that one is a linear 
function of the other. We have, in particular, analyzed the diagonal case where m% = m 2 , 
D m ,m, to be discussed below. For this family the chromatic polynomial P(D m m , q) has the 
form (2.9) with j max = 6 and a set of A's that are squares and cross products of those in 
(3.16). Since the equations defining the equality of magnitude of dominant A's can again be 
satisfied as 1/q — > 0, it again follows, by the criteria of [8j [20j [37] that the locus B extends 
infinitely far from the origin of the complex q plane. 

In a different direction, we have also studied a family F m (see below) for which the 
chromatic polynomial P(F m ,q) is of the form (3.16) with j max = 3, but with terms Xfj, 
j = 1, 2, 3, that are not simple polynomials, but instead are roots of a cubic equation, 
(18.16). In the m — > oo limit, the continuous accumulation set B for this family has bounded 
magnitude (does not extend infinitly far away from the origin in the q plane), as can be seen 
because the condition that defines B, namely the equality in magnitude of two dominant 
A^j's, cannot be satisfied for arbitrarily great |g|. 



IV. THE FAMILY B„ 



In Ref. PQ, we constructed and studied several one-parameter families of planar trian- 



gulations with chromatic polynomials of the multi-term form (3.16). (Some graphs in these 
families are listed in [27J-) in view of the general factorizations for the coefficients c G ( j 
that we have proved here, it is useful to express the results in terms of the reduced coefficients 
KG,j- For the family of bipyramid graphs B n , Eqs. (6.1)-(6.4) of Ref. [I] are equivalently 

= m + 2, which is well-defined for m > m m j 



3, with 



expressed as Eq. (3.16), with n(B m 
the reduced coefficients 

«B,1 = «B,2 = K B ,3 = 1 • (4.1 

V5 )/2 = 0.6180.. and, special case of (p.l8|, a# = 1 



1 



For this family, r(B 00 ) = 
Here and below we will index the members of the various families of planar triangulations as 
G p t,m, equivalent to the notation G pt>n that we used for some families in pQ via the relation 



V. THE FAMILY H m 

Here we remark on another family of planar triangulations H m p], which is well-defined 
for m > m m i n = 3 and has n(H m ) = m + 5. In Fig. [l]we show the lowest member of the 
family, H 3 . The next higher member, if 4 , is constructed by adding a vertex and associated 
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FIG. 1: Graph H 3 . 




FIG. 2: Graph H A . 

edges in the central diamond-like subgraph, as shown in Fig. [2j and so forth for higher 
members. We have given the chromatic polynomial for this family in Ref. pQ and have 



analyzed its properties there. In our present notation, P(H m ,q) has the form (3.16) with 

/^,i = (<?-3) 3 , (5.1) 
KH2 = q 3 - 9q 2 + 30g - 35 , (5.2) 



and 



Ktf,3 = -(?-3)(g-5) . 



(5.3) 
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FIG. 3: Graph L 4 . 



For this family, r = (7 - 3^ )/2 = 0.145898.. pQ. 



VI. THE FAMILY L r 



For comparative purposes, it is useful to study another family of planar triangulations 



with chromatic polynomials of the multi-term form (2.9). Here we denote this family as L m . 
It is well-defined for m > m min = 3 and has n(L m ) = m + 5. The lowest member of this 
family, L 3 , is the same as H 3 , shown in Fig. [TJ The next higher member, L 4 , is constructed 
by adding a vertex and associated edges in the central "diamond" , as shown in Fig. |3j and 
so forth for higher members. 

For this family of planar triangulations we calculate the chromatic polynomial P(L m , q) 



to be of the form (3.16) with 



k Lj1 = (g-2)(g-3) 2 , 
k l ,2 = q 3 ~ 9q 2 + 29q - 32 



(equal to \tc) an d 



k>l,3 = 2(g - 3) . 

P(L m , q) contains the factor P(K^ q) and has x{L m ) — 4. 
Evaluating P(L m , q) at q — r + 1, we find 



(6.1) 
(6.2) 

(6.3) 



P(L m , r + 1) = (-2 + 75 ) \{t - l) m + 2(r - 2) 



(6.4) 
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Consequently, 



r(L r 



-2 + 



1 - V5 



(6.5) 



As before, since |(1 — \fh )/2| < 1, the second term in Eq. (6.5) vanishes (exponentially fast) 
as m — > oo, so 

r(Loo) = -2 + V5 = 0.236068 (6.6) 

(to the indicated accuracy) and, as a special case of Eq. (3.18), a L = 1. 

We proved in general above that for any family of planar triangulations G pt}Tn with chro- 
matic polynomials P{G p ^ m) q) of the form (3.16), P{G p ^ mi q) has a zero that approaches 
t + 1 as m — )• oo. The families B m , H m , and L m (as well as others to be discussed below) 
illustrate this general result. In the present case, for odd m and hence even n = m + 5, this 
zero of P(L m , q) is slightly less than q — r+1, while for even m and hence odd n, the nearby 
zero is slightly greater than r + 1. We list these zeros in Table [T] for m from 4 to 15 (i.e., n 
from 9 to 20). 

If and only if m is odd, i.e., n is even, P(L m ,q) has another real zero somewhat larger 
than 3, which decreases monotonically toward 3 from above as m — > oo. As examples, for 
n = 8, 12, 16, 20, and 24, this zero occurs at approximately q = 3.61, 3.37, 3.25, 3.19, and 
3.16, respectively. 



VII. TWO-PARAMETER FAMILIES OF PLANAR TRIANGULATIONS, 



1)12 



In this section we introduce a substantial generalization to a two-parameter family 



Gpt, mi ,m2 °f planar triangulations involving the three A/s in (3.1), with a chromatic polyno- 
mial of the form 

3 3 

P(G ptjmi , m2 ,q) = ^2 C G,hi2 K 

Explicitly, 



V m 2 

\2 



(7.1) 



P(G pt> „, q) = c G , n (q ~ 2) mi+m2 + c G , 22 (g - 3) mi+m2 + Cg , 33 (-1 



\m\+m,2 



+ c G , 12 {q - 2) mi (<7 - 3) m2 + c G , 21 (g - 3) mi (<7 " 2) 



l«2 



+ c g , 13 (q - 2) mi (-l) m2 + c G , 3 i(-l) mi (g - 2) 
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TABLE I: Location of zero q z of P(L m ,q) closest to r + 1, as a function of the number of vertices, n = to + 5. 
Notation ae-v means a x 10"" here and in tables below. 



n 


Qz 


Qz - {T + 1J 


9 


2.630048 


0.01201 


10 


2.614750 


-0.003284 


11 


2.621594 


0.003560 


12 


2.616447 


-0.0015865 


13 


2.619250 


0.001216 


14 


2.617375 


-0.6588e-3 


15 


2.618477 


0.4432e-3 


16 


2.617774 


-2.598e-4 


17 


2.618200 


1.661e-4 


18 


2.6179335 


-1.005e-4 


19 


2.618097 


0.6294e-4 


20 


2.617995 


-3.8574e-5 



+ c G , 23 (g - 3) mi (-l) m2 + c G ,3 2 (-l) mi (<7 - 3) m2 . (7.2) 

(Below we shall often take i\ — i, ii — j to simplify the notation.) Clearly, if one keeps 
one of the indices m\ or m 2 fixed and varies the other, this defines an infinite set of one- 
parameter families of planar triangulations. For specific families G p t )mi ,m 2 we wm show how 



the general structure (7.1) reduces, in such cases, to the form (3.16) considered above for a 
class of one-parameter planar triangulations, with m being equal to the variable index, up 
to an appropriate integer shift. 

As before for the one-parameter families, an equivalent way to obtain the P{Gpt, mi ,m 2 , <l) 
is via a Taylor series expansion, in the auxiliary variables x\ and X2, of a generating function 
T(G pt , q, x\, X2). Equivalent to both of these is the property that -P(G ptimii . m2 , q) satisfies a 
two-dimensional recursion relation, for nil > (mi) min + 3 and m 2 > [jTi^) m in + 3, 

3 3 

P{G p t, mi ,m2y <?) + ^ ] ^ ] ^Gpt,hi2 P(Gpt,mi—ii,m2-i-ii = • (^-3) 
il=l 12 = 1 
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The coefficients in this recursion relation are given by 

3 3 3 3 



1 + e e ^4 2 = [ n ( x - a *i ^o] [ n ( x - ^ ■ ( ? - 4 ) 



H=l 12=1 ii=l 

Note that they satisfy the symmetry property 

We first derive a number of restrictions on the coefficients cq,^- As is true of any trian- 
gulation, P(G pt , mum2 ,q) = for q = 0, q = 1, and q = 2. The evaluation P(G pt , mu m 2 , 0) = 
reads 

c G ,n(-2) mi+m2 + c G , 22 (-3) mi+m2 + c G ,3 3 (-l) mi+m2 
+ c G , 12 (-2) mi (-3) m2 + c G , 21 (-3) mi (-2) m2 
+ CGil3 (-2) mi (-l) m2 + c G ,3i(-l) mi (-2) m2 

+ c G , 23 (-3) mi (-ir 2 + c G ,32(-l) mi (-3r 2 = at q = 0. (7.6) 

Since this equation applies for arbitrary m\ and m 2 in their respective ranges, it implies that 
c G,iii 2 = f° r an Hi ^2 at q = and hence that 

c G,iii 2 contains the factor q V ii, ? 2 . (7.7) 
It will often be convenient to extract this common factor, via the definition 

CGMi, = ^f- ■ (7-8) 
The evaluation P(G pttmi;m2 , 1) = reads 

[CG.11 + C Gi33 + C G ,13 + C Gl3 l](-l) mi+m2 + C G , 22 (-2) mi+m2 

+ [c G ,i2 + c G , 32 ](-l) mi (-2) m2 

+ [c G , 2 i + c G , 23 ](-2) mi (-l) m2 =0 at 9 = 1. (7.9) 
Since this equation applies for arbitrary m 1 and m 2 , it implies the conditions 

C G ,11 + C Gi33 + C Gj i 3 + C G)3 i = 0, C Gj22 = 0, 
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CG,12 + C G ,32 = 0, C G; 21 + c G23 = at q = 1 . (7.10) 

In particular, this implies that 

Cg,22 contains the factor q — 1 . (7-11) 
The evaluation -P(G p t, mi , m2 , 2) = reads 

C G ,llO mi+m2 + [c G , 22 + C Gi3 3 + C G)23 + C G , 32 ](-l) mi+m2 

+ [c G ,i2 + c G , 13 ]0 mi (-ir 2 + [c G , 2 i + c G ,3i](-l) mi m2 = at q = 2. (7.12) 

Since this equation applies for arbitrary mi and m 2 , including mi = m 2 = 0, it implies the 
conditions 

C G ,11 = 0, C G ,22 + C G ,33 + C G ,23 + C G , 3 2 = 0, 

c G ,i2 + c G ,i3 = 0, c Gj2 i + c G3 i = at q = 2 . (7.13) 

In particular, this implies that 

c G n contains the factor q — 2. (7-14) 

For families G pt ^ mitm , 2 with x(Gpt,m.i,m 2 ) — 4 for certain values of m x and m 2 , further condi- 
tions hold, as we shall discuss below. 

We next derive some further restrictions on the coefficients cq,^ from the condition that 
P(G p t jmi , m2 , q) must obey the Tutte upper bound when evaluated at q — r + 1. We consider 
families such that 

n(G ptjmi:Jn2 ) = mi + m 2 + (7.15) 
Carrying out this evaluation and calculating the ratio r(G pt , mi , m2 ), we have 



/ r _2\ m l+ m 2 / —1 \m 1 +m 2 

r(G pt , mum2 ) = (r - 1) -p c G ,n + c G , 22 (^ — yj + c G , 3 3^— J 

/r-2\ m 2 / 



/ -1 xm 2 

+ c G ,i 3 (— j +^il vT _ i/ 

(r - 2) mi (-1)™ 2 (-l) mi (r - 2) r 

+ C G,23^" TTZ-TZ" HC G , 32 - 



r — 


1 


T — 


2 


T — 


1 


-] 





(t - l)mi+m 2 ^ ( r _ !) 



mi+m 2 



(7.16) 
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The condition that r(G p t- mit m 2 ) — 1 f° r arbitrary mi and m2 implies that 

C G ,33 = C G ,13 = C G ,31 = C G ,23 = C Gj3 2 = at q = t + 1 . (7.17) 

By the same argument that we used above for the analysis of the coefficients of chromatic 
polynomials of one-parameter planar triangulation graphs, (7.17) implies that 

CG,hi 2 contains the factor q 2 — 3q + 1 if i\ = 3 or i 2 = 3 . (7-18) 

By taking either m 1 — >■ oo or m 2 —> oo, and requiring that the resultant ratio ^(G pii00)m2 ) or 
r (G P t, mi ,oo) must obey the Tutte upper bound, one deduces the inequality 

— l) 5 -^|c Gi n| < 1 atg = r + l. (7.19) 

We next generalize our result on a real chromatic zero that approaches q = r + 1 
for one-parameter planar triangulations to these two-parameter planar triangulations with 
P{Gpt^ mi ,m 2 i <?) °f the form (7.1 ). As in the p = 1 case, we assume that q is a real number in 
the interval [q w ,3). We will actually obtain two results, corresponding to m\ — > oo for fixed 
m2 and m 2 — > oo for fixed mi. For the first of these limits, the six terms proportional to 
cg,u, c g ,22, c Gj i2, c Gj2 i, c Gj i 3 , and c Gi23 all vanish (exponentially rapidly), so that 

P(G pt , mi ,m 2 ,q) ~ c G ,3 3 (-l) mi+m2 +c G ,3i(-l) mi (?-2) m2 +c G ,32(-l) mi (?-3) m2 for mi oo 

(7.20) 

But we have shown above in (7.18) that cc^^ contains the factor g 2 — 3g+l if = 3 or i 2 = 3. 
Since this factor vanishes at q = r + 1 in this interval [q w , 3), it follows that for sufficiently 
large mi, P(G p t >miim2 , g) has a real zero that approaches r + 1. With obvious changes, a 
corresponding argument shows that for sufficiently large m 2 and fixed mi, P{G p t. mi ,m 2 ,q) 
has a real zero that approaches r + 1. Clearly, the result also holds if both mi and m 2 get 
large. 



VIII. GENERAL FORM OF P{G pt ,, 



The generalization of our structural results for two-parameter families of planar trian- 
gulations, G pt , mi ,m 2 to j>parameter families is as follows. Let G vt ,m be a family of planar 
triangulation graphs involving the three A/s in (3.1) and depending on the p parameters 

i,i = l, -,P- Then 



m 



mi, 



,m r 



taking on integer values in the ranges m 8 > (m» 



P(G 



pt,m\,...m 



ii=i 
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In general, there are 3 P terms involving products of the A's (multipled by respective coeffi- 
cients) in this sum. 

This general form is of considerable interest. It shows that one can carry out a p-fold 
sequence of edge proliferations, each of which involves arbitrarily many additional edges, 
as indexed by the parameters mi,...,m p , with the chromatic polynomial P(G pij ^,g) still 
retaining the rather simple form (8.1) with the same set of three Aj's given in ( 3.1| . This 



is a much simpler situation than that in previous calculations of chromatic polynomials 
for multiparameter families of graphs. For example, in [28], Tsai and one of the present 
authors calculated the chromatic polynomial P(G ei ^ 2 ^ egim ,q) for a certain four-parameter 
family of cyclic chain graphs in which each subgraph on the chain has e\ edges above, and 
e 2 edges below, the main line, with e g edges between the subgraphs, and m subgraphs in 
all. Although the number N\ of Agj's for this family has the fixed value of 2, the A^j's 
have functional forms that depend on the parameters e%, e2, and e g (as is the case for the 
full Potts model partition function [29]). This property was also found to be true for (i) the 
chromatic polynomials P((Ch)k,m, C yc.,q) and P((Ch)k, m ,Mb.,Q) of cyclic (eye.) and Mobius 
(Mb.) strips depending on a homeomorphic expansion parameter k and the strip length, m, 
where N\ = 4 and three of the Xch/ S depended on k [20]; and (ii) P(H k>r ,q) for a family 
of "hammock" graphs H^^ with r "ropes" (linear sets of edges) joining two end vertices, 
with each rope having k "knots" (vertices), where again the N\ = 2 terms A#j depended 



on k and r [211 122]. The remarkable simplicity of the form (8.1 ) is a result of the restrictive 
property that G p t,m is a planar triangulation. We know that this simple behavior does not 
obtain even for the lowest case of one-parameter families for planar near-triangulations, from 
the explicit calculation the chromatic polynomials for free strips of the triangular lattice of 
length m and width L y = 2, 3 [31] (which are near-triangulations), where it was found 
that the A's changed with increasing width [31]. (Here a near-triangulation is defined as a 
graph such that all faces except one are triangles.) We also know that it does not hold for 
nonplanar triangulations, from explicit calculations of chromatic polynomials for the L y = 2 
[32] , L y = 3 [33], and L y = 4 [31] strips of the triangular lattice with doubly periodic 
(toroidal) boundary conditions. Thus, chromatic polynomials of multiparameter families of 
planar triangulation are especially amenable to exact analytic treatment. 

The P(G p t t m lt ...m p ,q) satisfy a p-dimensional recursion relation, for me > (me) min + 3, 
£ = 1, ...,p, namely 

3 3 

P(Gpt,m 1 ,...m v i q) + ^ ^ ' ' ' ^ ^GptM- ^p P{Gpt,mx—ii,...,m p —ipi Q) = 
ii=l ip=l 

(8.2) 
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where the b Gpuil ... ip are given by 



1 + e • • • e b ^~t* (ii ' - s n ri( i - a < *<) 



ii=l j D =l 



8=1 



1 i=l 



(8.3) 



Using the same methods as for p — 2, it is straightforward to generalize our results to this 
case, including (i) the conditions on the coefficients c Gpt i (where % = (ii...i p )) derived from 
the evaluations P(G p ^ m ,q) = for q = 0, 1, 2 and the Tutte upper bound at q — r + 1, 
and (ii) the results for r(G pttm ) and its limits as one or more of the m.j — > oo. Clearly, our 
result on a real zero in the interval [q w , 3) that approaches r + 1 also generalizes to this case 
of families G pt , m with p > 3. 

As in the p = 2 case, if one holds all but one of the mi, ...,m p fixed and allows one to vary, 



then the general form (8.1) reduces to (3.16) with m being equal to the variable parameter, 



up to an appropriate integer shift. 



IX. THE TWO-PARAMETER FAMILY D mum2 

We proceed to analyze our first explicit two-parameter family of planar triangulations, 
denoted D mxm2 (where D stands for the proliferation of a double set of edges). To explain 
the general method of construction of this family, we show in Fig. [4] the lowest member of 
the series, namely the graph D . We now add mi inner edges joining the uppermost vertex 
to the upper horizontal edge (thereby producing several such upper horizontal edges) and, 
separately, add m 2 inner edges joining the central vertex to the lower horizontal edge (thereby 
producing several such lower horizontal edges), with corresponding edges connecting to the 
lower central vertex. Thus, in the D miim2 graph, the uppermost vertex has degree 6 + mi, 
the central vertex has degree mi + m 2 + 4, and the lower central vertex has degree 4 + m 2 . 
To illustrate this, we show the graphs Di 2 , D 2 ^ in Figs. [5]and[6j 

From inspection of these graphs, it is evident how to construct -D mi m2 graphs with higher 
values of mi and m 2 . The number of vertices in the graph D mitm2 is 

n(D murrl2 ) = mx + m 2 + 9 . (9.1) 

For the chromatic number, we find that (i) if m 2 is odd, then x(D mi m2 ) = 4, and (ii) if 
m 2 is even, then x(D mi m2 ) = 3 if m x is even and x(-Dmi,m 2 ) = 4 if m x is odd. That is, 
denoting even as e and odd as o, x(D mijm2 ) = 3 for (mi,m 2 ) = (e, e) and x(D mijm2 ) = 4 
for (mi,m 2 ) = (e,o), (o, o), and (o, e). The chromatic polynomials for the D miyTn2 with 
X = 4 contain a factor P(K^, q) = q(q — l)(q — 2)(q — 3) (and some contain an additional 
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FIG. 4: Graph D 0fi . 




FIG. 5: Graph D h2 . 



factor of (q — 3)). Of course, all chromatic polynomials of triangulations have the factor 
P(K 3 ,q) = q(q-l)(q-2). 

By means of an iterative use of the deletion-contraction relation, we have calculated the 
chromatic polynomial P(D mitm2 ,q) for an arbitrary graph in this general two-parameter 
D mi ,m 2 family. We find that P(D mi ^ m2 ,q) has the form (7.1) with coefficients cjj^j that are 



rational functions of q. With the definition (7.8), we calculate 



Gnu 



(<Z-2) 



(9.2) 
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FIG. 6: Graph D 2 ,2- 



and 



(g - l)(g - 3) 5 (g 3 - 9q 2 + 30g - 35) 



C L>33 



q-2 

(g 5 - llg 4 + 46g 3 - 88g 2 + 74g - 23) (q 2 
(g-l)(g-2) 

CD,i2 = (g-2) 3 (g-3) 4 , 

c A 2i = (g-2)(g-3) 6 , 
(g-2) 3 (g 2 -3g + l) 



3g + l) 



CZ3,13 



CD,31 



CD, 23 



CD, 32 



g-1 

(g-2)(g 2 -3g 
g-1 

(<?-3) 4 (g-5)(g 2 - 



i: 



3g+l) 



g-2 

(g-3) 2 (g-5)(g 2 -3g+l) 



(9.3) 

(9.4) 

(9.5) 
(9.6) 

(9.7) 
(9.8) 
(9.9) 

(9.10) 



g-2 

It should be noted that the poles at q — 1 and g = 2 in certain of these Cp^- coefficients 
are cancelled in the actual evaluation of P(D mi>m2 , g), which is, as it must be, a polynomial 
in g. Furthermore, not only are these poles cancelled, but also the resultant P(D mi)7n2 ,q) 
vanishes at g = 0, g = 1, and g = 2. It is easily verified that the coefficients co,ij satisfy the 
requisite conditions (7.10) and (7.13) for these zeros to occur. 
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Furthermore, since cd, 11 = 0^,33 = cd,iz = cd,3i = 3/2 with the other cn,ij = at q = 3, 
it follows that 

P(D mi>m2 ,3) = l[l + (-l) mi+m2 + (-1)™ 2 + . (9.11) 

This vanishes for (mi, m 2 ) = (e, o), (o, o), and (o, e), and is nonvanishing for (mi, m 2 ) = (e, e) 
in our notation above, in agreement with our result on the chromatic number x(D mi m2 ). 

We proceed to discuss the evaluation of the chromatic polynomial P(D mirn2 , q) at q = r+1 
and the comparison with the Tutte upper bound. We have 



P(D mum2 ,T + 1) = (9 - 4^5 )(t - I)"'*"* + ( 445 2 199 ^ )(t - 2)"«- 

-38 + 17V5N (T _ ir(T _ 2r 



-rrt 2 



+ 



+ 



2 / 
•199 + 89^5 



)(r-2) mi (r- l) ma (9.12) 



2 7 

Comparing this with the Tutte upper bound (r — 1 ) m i+ m 2+4^ we ij ave 

3-V5 ^65 - 29^ \ /l - v 7 ^ \ mi+m2 
^-11 + 5^ \ ^1 - \ m2 



+ ,-29+13^ )( l-^ r (913) 



We list a number of values of r(D mitm2 ) in Table |H| One may investigate the behavior of 
^(Dm 1 ,m 2 ) as mi —> 00 for fixed m 2 and as m 2 — > 00 for fixed mi. Because the quantity 
(1 — )/2 that is raised to the powers indicated in (9.13) is negative, it follows that, if one 
keeps m 2 fixed and increases mi, then r(D mim2 ) does not approach r(D oo m2 ) monotonically, 
although the members of the subsequences r(D mi<m . 2 ) with even (odd) m 1 approach r(D oo jn2 ) 
monotonically from above (below), respectively. Similarly, if one keeps mi fixed and increases 
m 2 , then r(D mitrri2 ) does not approach r(D mitOQ ) monotonically, although the members of the 
subsequences r(D mijm2 ) with even (odd) m 2 approach r{D mi)OQ ) monotonically from above 
(below), respectively. The results are 

,^ /„ 3-V5 /-ll + 5^ \ /l - V5 \ m 2 
rpoo, m2 ) s Jim o r( J D mi , m2 ) = + ( — V — J (— ^~ ) (9-14) 

KAni.oo) = ^^/(Arn.mJ = g + V 2 ) \ 2 / ^ ^ 
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In the limit where one takes both mi and m.2 to oo, one has 



r(D c 



lim lim r(D mim2 ) = lim lim r(D mit 



mi->oo m2->oo 



mi— >oo mi->oo 



3 - V5 



2- T = 0.381966... 



(9.16) 



As 77i2 increases from to oo, r(D oom3 ) decreases (non-monotonically) from the value 

r(£>oo,o) = -4 + 2V5 = 0.4721359... (9.17) 



to the value in (9.16), and as m x increases from to 00, r(Z? mi)00 ) decreases (non- 

-13 + 6^ = 0.4164078... 



monotonically) from the value 



(9.18) 



to the value in (9.16). As a consequence of the relation (9.34) (see below), it follows that 

r{D 00tk+2 ) = r{D Koo ) (9.19) 

In general, the maximal value of r{D mi)Tn2 ) occurs for the member of the D mi)7n2 family with 
the minimal values of mi and 7772, namely for Do,o- This property is similar to the property 
that the maximum value of r(G pt ) for all planar triangulations G pt occurs for the G pt with 
the minimum number of vertices, namely the single triangle, K 3 . 



We next show that our general form for P(D miiTn2 , q) reduces to (3.16) when either 7772 is 



XT 



(9.20) 



held fixed and mi varies, or vice versa. If we keep 777,2 fixed and vary mi, then we can write 

3 3 

P(D mum2 ,q) = J2 [J2 C ^ X T 
i=i 3=1 

The sum Y^j=i c D,ij XJ 12 contains a factor Xj, which we combine with the A™ 1 , to make A™, 
where 

m = m 1 +A. (9.21) 



This shows that P(D mi;Tn2 ,q) has the form (3.16) with m given by (9.21); explicitly, 



P{D m -A,m 2 ,q) - J^CD mvU) ,i A 



in 

i 1 



i=l 



where 



CD 



3=1 



771.2 

3 



(9.22) 



(9.23) 
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TABLE II: Values of the ratio r(D m i m 2). The rows and columns list mi and mi, respectively, so that, for 
example, r(D lj2 ) is the entry 0.3769. 



mi, m-2 


u 


i 
i 


o 
Z 


Q 
O 


A 
1 








oo 





0.5836 


0.3131 


0.4803 


0.3769 


0.4408 


0.4013 


0.4257 


0.4164 


1 


0.4033 


0.3344 


0.3 my 


0.3506 


0.36by 


0.35oo 


0.3631 


U.360 ( 


2 


0.5147 


0.3212 


0.4408 


0.3669 


0.4126 


0.3843 


0.4018 


0.3951 


3 


0.4458 


0.3293 


0.4013 


0.3568 


0.3843 


0.3673 


0.3778 


0.3738 


4 


0.4884 


0.3243 


0.4257 


0.36305 


0.4018 


0.3778 


0.3926 


0.3870 


5 


0.4621 


0.3274 


0.4106 


0.3592 


0.3910 


0.3714 


0.3835 


0.3789 


6 


0.4783 


0.3255 


0.4200 


0.3616 


0.3977 


0.3754 


0.38915 


0.3839 


oo 


0.4721 


0.3262 


0.4164 


0.3607 


0.3951 


0.3738 


0.3870 


0.3820 



Here, since these coefficients depend only on m 2 , and not on mi, we have introduced the 
notation D m2 ig\ to refer to the entire family D miJn2 with fixed m% and variable mi, where I 
indicates that m 2 describes the edge proliferation in the lower part of the graph. Expressing 
(3.13 )-(3.15 ) in our notation (and suppressing the q arguments), we have 



c Dm2(l) ,2 = q{q - 1) Kb, 



and 



CD 



2 CO' 



q(q 2 - 3q + 1) K Dr> 



2(<0> 



Similarly, if we keep mi fixed and vary wi2, then we can write P(_D mijm2 , q) as 



3 3 
j=l i=l 



(9.24) 
(9.25) 

(9.26) 
(9.27) 



The sum Yli=i c D,ij A™ -1 contains a factor A|, which we combine with the A™ 2 , to make A 
where, for this one-parameter reduction, 



in 

j ' 



m = m 2 + 2 . 



(9.28) 
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This shows that P(D mi ^ m2 ,q) has the form (3.16) with m given by (9.28); explicitly, 



where 



P(D mi , m -2,q) - 22 c D mi(u) ,j K 

3=1 



\— 2 \ \mi 

i=l 



(9.29) 



(9.30) 



Here again, since these coefficients depend only on mi, and not on m 2 , we have introduced 
the notation D mi ( u ) to refer to the entire family D mi>m2 with fixed ni\ and variable m 2 , 
where u indicates that mi describes the edge proliferation in the upper part of the graph. 
As before, we write 

C D mi(u) ,l = ?KD mi(u) ,l , (9-31) 

°D mi(u) ,2 = q(q ~ 1) Ka» i(i0 ,2 > ( 9 - 32 ) 



and 



Cd , , 3 = q(q —3q + l)KD M 3- 

^m^u)! -iV^ ^ / ^ mi(u) i° 



We find that 



and thus 



CD. 



fe+2 



CD fe(u) ,i for z = 1, 2, 3 



However, we note that for arbitrary q, 



for i = 1, 2, 3 . 



(9.33) 
(9.34) 
(9.35) 



P(D mi>m2 ,q) ^ P(D m2tTni ,q) unless mi = m 2 



(9.36) 



X. THE FAMILY D m _ 4 o 



We proceed to examine a number of different D muJn2 families of planar triangulations, 
with m 2 held fixed. Then, we will analyze analogous families with mi held fixed, and finally, 
we will investigate families in which both mi and m 2 vary together, and are related in a linear 
manner. For a given graph D mitm2 , one can use either our general result for P(D mi>m2 ,q) 



above or either of the one-parameter reductions, (9.22) or (9.29). However, we shall be 



interested in the limits mi — > oo with m 2 fixed, and m 2 — > oo with mi fixed, and, to study 
these, it is convenient to use the one-parameter reductions of our general formula. 
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TABLE III: Values of the ratios r(D mit00 ) and r(D 00jjn2 ). Note that r(Dk :00 ) = r(D OQ: k-2) for k > 2. 



r( 7") \ r( D \ 
i \ J -'oc,m2)i 1 y^mi ,oo ) 


clilcliy L1C 


IlUlIlcl lCdl 


^(^>oo,o) 


-4 + 2y/5 


0.472136 


K^oo.l) 


f -1 r* 1 r—7 ft* \ 1 c\ 

(-15 + 7V5 )/2 


0.326238 


r{D 00,2) = r(iJo,oo) 


— 13 + 6\/5 


0.416408 


r(I?oo,3) = r{Di }0O ) 


-22 + 10\/5 


0.360680 


r(-Doo,4) = f*(L>2 i00 ) 


(-73 + 33V5 )/2 


0.395122 




-60 + 27a/5 


0.373835 


^(,^00,6) — ^(^4,ooj 


no 1 a a fri 
— \)o + 44 v 


u.oooyyi 




(-319 + 143^/5 )/2 


0.378860 


KAjo.s) = r(L> 6j00 ) 


-259 + 116^/5 


0.383885 


KA>o,9) = r ( D 7,oo) 


-420 + 188^5 


0.380780 


KAx),lo) = ?"(£>8,oo) 


(-1361 + 609^5 )/2 


0.382700 


^(-^00,00) 


(3- V5)/2 


0.381966 



We begin with a study of the family D mi:0 = D m _4 i0 with mi > 0, i.e., m > 4. From 
(9.1), we have n(_D m _ 4) o) = m + 5. For the coefficients that enter into the equation (3.16), 



our general formulas (9.23)-(9.26) yield 



n Do{e)i3 = q 3 -9q 2 + 29q-32 



(10.1) 



(equal to Xtc)- Because these coefficients K>D 0(t) ,j are all the same, Xtc is a common factor, 
so for all m, the three zeros of Xtc are zeros of P(_D m _4 5 o, Of these, one is real, namely 



q w , given in (3.25). In accordance with our general analysis above, P(D Tn _4 j0 ,q) also has 
a zero, denoted q z , that approaches r + 1 as m increases. We list this zero for m — 1 to 
m = 16 in Table IV As is evident from this table, for odd (even) m, this zero is slightly 



above (below) r + 1. 
For the evaluation at q 



t + 1, we have 



P(An-4,o,r + l) = (-4 + 2^ )(r-l) m + (3- V5 )(r-2) 



(10.2) 
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TABLE IV: Location of zero q z of P(D m _4fi,q) closest to r + 1, as a function of the number of vertices, 
n = m + 5. 



Tl 


yz 


— (t 4- L) 


1U 


Z.b I /olo 


u.uoy < o 


11 


Z.oy4ozy 


— U.UZoZl 


iz 


z.bobllo 


U.UloUo 


lo 


z.bUyioU 


— U.oyU4e-/ 


1 A 
14 


z.bZ4oob 


u.bozze-z 


10 


O fil A KA 1 
Z.bl4041 


— u.o4yoe-z 


lb 


z.bzUoob 


Z.oZZe-o 


1 7 
1 / 


O £1 ££VQ 

z.blbb i o 


— l.oble-o 


lo 


z.bloyUo 


U.o < loe-o 


19 


2 61 7509 


—0 5254p-3 


20 


2.618364 


0.3301e-3 


21 


2.617832 


-2.017e-4 


22 


2.618160 


1.2560e-4 


23 


2.617957 


-0.7725e-4 


24 


2.618082 


0.4790e-4 


25 


2.618004 


-2.954e-5 



so that 

1 /5 

r(D m . Afi ) = -4 + 2v/5 + (3 - ) ) 

Hence, 

r(D OOi0 ) = -4 + 2^5 = 0.4721359... 

with a DQii) = 1. 



(10.3) 
(10.4) 



XI. THE FAMILY D 



m— 4,1 



We continue with a study of the family D 



,„,.! = D m _4,i. From plj ), we have n(D m _ 4A ] 

as 



m + 6. Our general formulas (9.23 )-(9.26 ) give the coefficients KDunj 

KD 1(t) ,i = (q~ 3)(? 3 - 9g 2 + 30g - 35) , 



kd 1( «,2 = g 4 - 12g 3 + 58g 2 - 133g + 119 



11.1) 



11.2) 



and 



KD m , 3 = -{q - 3)(2g 2 - 14g + 25) . (11.3) 

If m is even, then P(_D m _ 4j i, g) has not only the factor q(q — l)(q — 2)(q — 3), but also an 
additional factor of (q— 3). According to our general analysis above, P(D m _4 i i, q) has a real 
zero that approaches r + 1 as m — > oo. We also derived the result that for sufficiently large 



m, a chromatic polynomial of the form (3.16) has another real zero in the interval [q w ,3) if 



and only if KG pt ,3 has a zero in this interval. For the present family, kd 1w ,3 has zeros at q = 3 
and the complex-conjugate pair q = (7±z)/2, but does not have a zero in the interval [q w , 3), 
in accordance with the fact that P(Z) m _4 i i, g) also does not have a zero in this interval. 
For the evaluation at r + 1, we compute 



P(D 



m— 4,1) 



r + 1) 



25 + 11^ 



l) m + (-9 + 4v / 5)(r-2) T 



so that 



and 



r(L> 



m-4,1; 



-15 + 7^ /ll-5 v / 5wl 
2 h V 2 J 



x/5 



r(D, 



■15 + 7^5 



00,1/ 



0.3226238 



11.4) 



11.5) 



11.6) 



with a D = 1. 



XII. THE FAMILY D 



m-4,2 



We next study the family £) mi ,2 = -Dm-4,2 with mi > 0, i.e., m > 4. Note that is the 
same as the graph denoted G ce i2 in Fig. 8 of [T]. From (9.1), we have ra(D m _ 42 ) = m + 7. 
For this family our general results give 



K Dm> i = q- 15g 4 + 94g 3 - 303g 2 + 498g - 332 



K D 2(e) ,2 



q 5 - 15g 4 + 95g 3 - 317g 2 + 553g - 39£ 



12.1) 
12.2) 
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and 

k d 2W ,3 = ~{q 4 ~ 16g 3 + 91g 2 - 225g + 206) . (12.3) 

In Table [V] we list (real) zeros of P(D m _4 ! 2,q) in the interval q G [q w ,3) as a function of 
n. As proved above, one zero approaches r + 1 as m — > oo. In the same limit, our general 
analysis above shows that P(D m _ i 2 ,q) has real zero(s) in the interval [q w ,3) correponding 
to the zeros of «d 2 ,«,3 hi this interval. This quartic polynomial has a zero at 

q = 2.7227000945... (12.4) 

together with one more real zero at q — 6.955106.., outside the interval [q w , 3), and a complex- 
conjugate pair. Hence, P(D m _ 4j2 , q) has another zero in the interval [q w , 3), which is present 



for m > 5, and this approaches the zero of kd 2W ,3 given in (12.4) as m — > oo. For even 
m > 6, i.e., odd n > 13, the real zero near to this asymptotic value (12.4) increases toward 
it, while for odd m > 5, i.e., even n > 12, the nearby real decreases toward the asymptotic 
value. As noted above, the graph D 2 coincides with the graph Gqm,i of PQ, for which there 
is no real zero close to r + 1; instead, the zeros closest to r + 1 comprise a complex-conjugate 
pair at q = 2.641998 ±0.014795z. Correspondingly, for D 2 there is no zero q' z in the interval 
[q w ,3). For all of the D m -4,2 with m > 5 in Table [VJ the first real zero q z in the interval 
[q w , 3) is, in fact, the closest to r + 1. 
The evaluation at q = r + 1 yields 

-69 + 31^ 



so that 



Hence, 



P(D m ^, 2 ,T + 1) = ( i— ^- ) (t - 1)"' + (27 - 12V5 )(t - 2)'» , (12.5) 

r(A„-«) = -13 + 6^5+ (Hi^5 , y _ 



r(D 0O)2 ) = -13 + 6^ = 0.41640786... (12.7) 



with a Dm 



XIII. THE FAMILY L> m _ 4i3 
The final family that we study in this series is D mi 3 = D m _4 3 , with m\ > 0, i.e., m > 4. 



From (9.1), it follows that the graph D m _ i:3 has n(D m _4 5 3) = m + 8. The graphs and 
I?2,3 are shown in Figs. [7] and [8} 
For the coefficients ^D 3(t) ,j we have 

KD m ,i = {q- 3)(g 2 - 5g + 7)(g 3 - 10g 2 + 38g - 49) , (13.1) 
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TABLE V: Location of real zeros of P(D m _4 2, q) in the interval q £ [q w , 3), as a function of the number of 
vertices, n = m + 7. Here the notation nz means that there is no second real zero in the interval [^,3). 



Ti 






11 


c.c. pair 




nz 


12 


2.614614 


2 


818897 


13 


2.621801 


2 


689610 


14 


2.616506 


2 


762806 


15 


2.619226 


2 


705035 


16 


2.6174035 


2 


741044 


17 


2.618462 


2 


713055 




Z.Ol (IOU 


2 


731543 


19 


2.618194 


2 


717464 


20 


2.617938 


2 


727100 


21 


2.618094 


2 


719886 


22 


2.617997 


2 


724931 


23 


2.618057 


2 


721202 


24 


2.618020 


2 


723845 



K Dma = g 6 - 18g 5 + 141g 4 - 613g 3 + 1551g 2 - 2152g + 1271 , (13.2) 

and 

«jj 3(<)l3 = -(? - 3) 2 (g 3 - 12g 2 + 48g - 67) . (13.3) 



In Table VI we list real zeros of P(D m _4 i 3,g) in the interval q G [q w ,3). As is evident 
in this table, in addition to the zero that approaches r + 1, there is a second real zero in 
this interval if and only if m (and hence n) is even. We can prove that, for the subset of 
D m _4 )3 with even m where this second zero is present, it approaches q = 3 from below as 
m — > oo. The proof is as follows. According to our little theorem above on real chromatic 
zeros besides the one (or complex pair) that approach r + 1 as m — > oo, there is a second 
zero in the interval q G [q w , 3) if and only if k g ^ has a real zero in this interval. Now for 
the D m _4 i 3 family, ft£> 3( «,3 has no real zero in the interval [q w , 3). (Its zeros are at q = 3, with 
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FIG. 7: Graph D li3 . 




FIG. 8: Graph D 2fi . 

multiplicity 2, at q = 5.44224957... and at q ~ 3. 278875± 1.249025.) Hence, according to our 
little theorem, as m — » oo on even integers, in addition to the real zero that is near to r + 1, 
the other zero must approach 3, so that in this limit, there is no other zero in the interval 
[q w , 3). This family may thus be added to the two known families given (in his Theorem 4) 
by Woodall in Ref. p2] as examples of one-parameter families of graphs, each of which has 
a chromatic zero that approaches 3 from below as the parameter (m here) goes to infinity. 
The evaluation at q = r + 1 yields 

P(An-4,3, r + 1) = (94 - )(r - l) m + (-76 + 34^ )(r - 2) m , (13.4) 
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TABLE VI: Location of real zeros of P(D m _4^, q) in the interval q G [q w , 3), as a function of the number 
of vertices, n = m + 8. Notation nz means that there is no second zero in this interval. 



Tl 






12 


2.614614 


2.818897 


13 


2.619530 


nz 


14 


2.616973 


2.847527 


15 


2.618625 


nz 


16 


2.617649 


2.866268 


17 


2.618264 


nz 


18 


n />i to on 

2.617889 


2.880165 


19 


2.618122 


nz 


on 


9 R1 7Q7Q 


z.oyuyoo 


21 


2.618068 


nz 


22 


2.618013 


2.899700 


23 


2.618057 


nz 


24 


2.618020 


2.906905 


25 


2.618039 


nz 


26 


2.618031 


2.912980 



so that 

r(D m _ 4 , 3 ) = -22 + 10^5 + (18 - 8^ ) ) . (13.5) 

Hence, 

r(£>oo,3) = -22 + lOv^ = 0.36067977... (13.6) 

with a D = I. 



XIV. THE FAMILY L> , m _ 2 
As an example of families with mi fixed and variable m 2 we discuss the family -Do, m2 



A),m-2- A graph in this family has n(D m _ 2 ) = m + 7. In accord with our result (9.35) 



KD 0(u) ,i = KD 2(u) ,i for i = 1, 2, 3 . (14.1) 

P(Do, m -2, <?) has a real zero near to r + 1, which approaches this point as m — > oo. Fur- 
thermore, since the coefficient «d , ,3 has a real zero in the interval [q w ,3), at the value in 



(12.4), it follows from our general analysis above that for sufficiently large m, P(D m _ 2 , q) 
has a real zero that approaches this value. These zeros approach their respective values in a 
manner similar to that discussed for the family -D m _ 42 . 



XV. A SYMMETRIC TWO-PARAMETER FAMILY S mitTna 

In this section we study a two-parameter family of planar triangulations S miiTn2 which are 
symmetric under interchange of the parameters: 

&m\,m 2 S m2 ^ mi . (15.1) 

We show the lowest member of this family, S^q in Fig. [9] and another member, 2 — £2,1 i n 



Fig. 10 From these it is clear how to construct the general graph S mitrri2 in this family. We 
have n(S muma ) = mi + m 2 + 7. The chromatic number is x(S mit m 2 ) — 4, and P(S mi ^ m2 ,q) 
contains the factor P(K^ q) = q(q — l)(q — 2)(q — 3). 



We have calculated P(S mi)m2 , q) and find that it involves the same three A's as in (3.1) 



but with an interestingly different form than P(D mi:ni2 , q). Given the symmetry (15.1), it 
follows that the coefficients in P(S mi , m2 ,q) satisfy 

cs,ij = c s ,ji • (15.2) 

Consequently, although there are nine terms of the form A^AJ 12 in P(S mi ^ m2 ,q), there are 
only six independent coefficients cs,ij to begin with, and we find that two of these vanish, so 
that there are only four independent, nonvanishing coefficients cs^j- Explicitly, we calculate 

cs,ij = qc s ,ij , (15.3) 

with 

CS,U = C S ,13 = C S ,31 = , (15.4) 
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FIG. 9: Graph S ,o- 




CS,22 — 



0?-i)0?-3) f 



and 



9-2 

C5,i2 = c 5 ,2i = (g-2) 3 (g-3) 2 , 
(g-3) 2 (g 2 -3g+l) 

CS,23 = Cs,32 = ~ 

(g-l)(g-3)(g 2 -3g + l) 



C5,33 



9-2 



(15.5) 
(15.6) 
(15.7) 

(15.8) 
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so that 



_L r- < \ m i \ m 2 I \ rrn \ m 2 \ , „ Ami \rri2 1 \mi \m2^ 
+ Cs,12{^i A 2 + A 2 A 1 ) + CS,23{A 2 A 3 + A 3 A 2 



(15.9) 



As before, the poles cancel in the calculation of P(S mitm2 ,q) and, furthermore, these co- 
efficients satisfy the requisite identities so that P(S mi ^ m2 ,q) = for q = 1, 2, 3. These 
are special cases of (7.10) and (7.13) that incorporate the properties that cs t ij = csji and 



e<?,n = c 5i i 3 = c s>31 = 0, namely, 

CS,22 = 0, C5-33 = 0, 05,12 + C S ,23 = at q = 1 

and 

cs,22 + 2c Sj i2 + c 5i3 3 = 0, c s ,i2 = at q = 2 . 



As a special case of (7.17) we also have 



CS,23 = cs )32 = c 5i33 = at q = t + 1 . 
Finally, the condition P(S mi)m2 ,3) = is equivalent to 

cs,33 = at q = 3 



(15.10) 
(15.11) 

(15.12) 
(15.13) 



For the comparison of P(S mi>m21 q) at q = r+1 with the Tutte upper bound (r— 1 ) m i+ m 2+2^ 
we have 

1-^/5 ^«ii+m 2 



r(£ TO1 , m2 ) = |(9-4^)( o 

1 - \ mi /l - V5 



+ (-2 + V5 



(15.14) 



This decreases (non-monotically) in magnitude as mi increases for fixed ra 2 and as m 2 
increases for fixed m 1; approaching zero exponentially rapidly as either of these parameters 
goes to infinity. Thus, 



lim r(S mi 



mi->oo 



TT12I 



lim r(S mi 



m2->oo 



m.2 j 



. 



;i5.i5) 



It is also of interest to analyze the one-parameter reductions of P(S mitm2 ,q) for variable 
mi and fixed m 2 and vice versa. These yield identical results, because of the symmetry 



(15.1). Hence, without loss of generality we consider variable m 1 and fixed m 2 and find that 
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in this case P(S mijTn2 , q) reduces to (3.16) with m = mi + 2, which we write as P(S m -2,m 2 , <?)• 
Since the coefficients only depend on m 2 and not m, we denote them by cs m2 ,i- They are 
given by 



Thus, in terms of the corresponding k,s„ 



A7 



lll 2 



K Sthl = (q-2)(q-3) 2 , 

k So ,2 = Xtc = q 3 - 9g 2 + 29q - 32 
^o,3 = 2(q - 3) . 



(15.16) 

(15.17) 

(15.18) 
(15.19) 



and so forth for higher values of m 2 - 

This family exhibits a number of interesting properties. Among these is the fact that 
out of the possible 3 2 terms in (8.1) for p = 2, some may be absent because of vanishing 
coefficients cc,ij- In particular, the term c^uA™ 1 2 that would be dominant in the limit 
where the parameters mi — > oo and m 2 — > oo, may be absent, so that in this limit, r(G 00 ^ 00 ) 
may be zero. 



XVI. FAMILIES OF THE FORM G 



pt,m\ ,rri2 



WITH mi = m 2 



In previous sections we have analyzed the chromatic polynomials of special cases of two- 
parameter families of planar triangulations G p t, mi , m2 as a function of mi with m 2 held fixed, 
and vice versa and shown how they reduce to (2.9) with j max = 3. A different type of 
special case in which G pt) m ljrn2 reduces to a one-parameter family is obtained by requiring 
that mi and m 2 be linearly related to each other. The simplest such example of this type 
of reduction is the diagonal case obtained by requiring that mi = m 2 . For general families 
G 



pt,m\ ,m,2 



that satisfy (7.15) and for which P{G 



pt,77X\ ,7772 1 



n(Gpt A k) = 2k + (3 and that P{G 



pt, m\,m\ i 



, q) is of the form (7.1), it follows that 

= 6. We 



q) reduces, to the form (2.9) with j 



use the shorthand Gd to denote a generic G pt , 



We have 



P(G d<mumi ,q) - y^c Gd ,j(A GdJ ) r 
where m = mi + 5m, with 5m depending on the family, and 



(16.1) 



\\ = {q-2f 



(16.2) 
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\ Gd ,2 = \l = (q - 3) 2 , 

Ac d ,3 = A3 = 1 , 

AG d ,4 = A 1 A 2 = (g-2)(g-3) 
A Gd ,5 = A1A3 = -(g-2) , 



and 



A Gd 



: A 2 A 3 = -{q - 3) . 
The corresponding coefficients are 

CG d ,j = qc Gd ,j forj = l,..,6, 



with 



and 



CG d ,l = cg,u , 
CG d ,2 = Cq,22 , 
CG d ,3 = Cg,33 , 
CG d ,4 = C G ,12 + CG,21 , 
CG d ,5 — CG,13 + Cg,31 , 

CG d ,6 = Cg,23 + Cg,32 ■ 



(16.3) 
(16.4) 
(16.5) 
(16.6) 

(16.7) 

(16.8) 

(16.9) 
(16.10) 
(16.11) 
(16.12) 
(16.13) 



(16.14) 

The coefficients cc d ,i, % = 1,...,6 satisfy various conditions that follow from those that we 



have derived for the coefficients CQ,ij in (7.7), (7.10), (7.13), and (7.17). These are 

c Gd ,i + c Gd ,s + c Gd , 5 = 0, C Gd ,2 = 0, 

cc d ,A + c Gd ,6 = at q = 1 , 

C Gd ,l = 0, C Gd ,2 + C Gd ,3 + C Gd ,6 = 0, 



and 



Hence, 



and 



cc d A + c Gd ,6 = at q = 2 , 
CG d ,3 = c Gd , 5 = c Gdfi = at g = T+l. 

°G d ,2 contains the factor q — 1 , 
c G d ,i contains the factor q — 2 , 

CG d ,i contains the factor q 2 — 3q + 1 if i = 3, 5, 6 



(16.15) 

(16.16) 
(16.17) 

(16.18) 
(16.19) 

(16.20) 
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XVII. THE FAMILIES D mi ,m 2 AND S mi ,m 2 WITH mi = m 2 



We now discuss two explicit examples of the diagonal special case of a two-parameter 
planar triangulation family, namely D mum2 and S mi>m2 with m\ = m 2 . We shall use the 



shorthand notation and Sd to refer to these entire respective families. From (9.1), we 
have n(D mumi ) = 2m x + 9 and n{S mumi ) = 2m\ + 7. 



The chromatic polynomial P(D mi>mi ,q) has the form (16.1) with j 



6 and 5m = 0, 



i.e., m = mi. The lowest member of this family, the graph -Do,o> was shown in Fig. [4] and 
Z^2,2 was shown in Fig. [6j The chromatic number x(D m ,m) is 3 if m is even and 4 if m is 
odd. The coefficients are 

CD d ,i = c D)ii for 2 = 1, 2, 3, 

CD d ,4 = ci 2 + c 21 = (q - 2){q - 3) 4 (2g 2 - lOq + 13) , 
(q - 2)(q 2 - 4g + 5)(g 2 - 3g + 1) 



CD d ,5 = C13 + C31 



and 



c D d ,6 — c 23 + c 32 



9-1 

(g - 3) 2 (g - 5)(g 2 - 6g + 10)(g 2 - 3g + 1) 
q-2 



(17.1) 
(17.2) 

(17.3) 



Since x(D m ,m) = 4 for odd m, it follows that 

c Dd ,i ~ c Dd:3 + c Ddi5 = for q = 3 . 

We calculate 

'47-21^5 



(17.4) 



;i7.5) 



P(D mjm ,r+l) = (r + 1) 



r - 1 



2 m 



+ 



/1165 - 521^5 



) ( T _ 2 ) 2m + (-360 + 161 V5 )[(r - l)(r - 2)]' 



(17.6) 



The ratio r(D m ^ m ) is 

3-^5 /(>•'> 29 v / ">\/l v'-"' ^ 



+ (-20 + 9^ 



1 - V5 



Hence, defining r(Dd, oo) = linim^oo r{D m ^ m ), we have 

r(D d , oo) = 3 ~ — = 0.381966. 



(17.7) 



(17. 
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and dD d = 1- 

In contrast, the chromatic polynomial P(S miimi ,q) has the form (2.9) with j max = 4. 
The lowest member of this family, the graph S^o, was shown in Fig. [9j The chromatic 
number x(S m ,m) = 4. The expression for P(S m)mi q) follows immediately from (15.9) and 

has jmax 4, 



P(S m>m , q) = c 5 , 22 (A 2 ) im + c 5 ,33 + 2c s ,i 2 (A 1 A 2 ) m + 2c s , 2 3(A 2 A 3 )' 



(17.9) 



where we used the fact that (A; 
result that r(S , 00j00 ) = 0. 



\2m 



1. The ratio r(S m ^ m ) follows from (15.14), with the 



XVIII. THE FAMILY F„ 



In this section we construct and study a family of planar triangulations, denoted F m , 

= 3, but the Xpj are not 



with the property that P(F m ,q) has the form (2.9) with j 



given by (3.1), but instead are roots of a certain cubic equation. The number of vertices is 
n(F m ) = m + 4. This family is useful as a contrast to the other one-parameter families of 



planar triangulations with chromatic polynomials of the form (2.9) that we have constructed 



The construction of members of this family is somewhat more complicated than that of the 
other families analyzed in this paper, and accordingly, for illustration we include several 
graphs, namely F m with m — 3, 4, 5, are shown in Figs. [Tl|fT3| As these show, starting 



from a given member F m , one constructs the next higher member F m+1 in an interleaved 
manner, first adding a new set of edges one of which emanates from the lower left-hand 
vertex of the graph, and then a new set of edges one of which emanates from the uppermost 
vertex, and so forth. We note that in contrast to the previous planar triangulations with 



chromatic polynomials of the form (3.16), the degrees of the vertices remain bounded as 
m — > oo for this family. 

As in earlier works p-[TU]. it is most convenient to express the P(F m , q) via a generating 
function, T(F, q, x), which is a rational function in q and an auxiliary expansion variable x, 
of the form 

Af(F,q,x) 



T(F,q,x) 



where the numerator and denominator are 



V(F,q,x) 



Af(F, q, x) = a F Q + a F iX + ap j2 



x 



and 



V(F, q, x) = 1 + b FA x + b F)2 x 2 + 6 Fj3 x 3 



18.1 



'18.2) 



18.3) 
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FIG. 11: Graph F 3 . 




FIG. 12: Graph F 4 . 

with dF,j and &fj being polynomials in g. The chromatic polynomial P(F m ,q) is then given 
as the coefficient in the Taylor series expansion of this generating function: 

oo 

r(F,q,x) = J2P(F m+ i,q)x m (18.4) 

m=0 

Using an iterative deletion-contraction method, we have determined this generating func- 
tion. We find 

a F , = q(q-l)(q-2)(q-3) 2 (18.5) 
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FIG. 13: Graph F 5 . 




FIG. 14: Graph F 6 . 

a Ftl = q(q-l)(q-2)(2q-5) (18.6) 

a F , 2 = q(q-l)(q-2) 2 (q-3) 2 (18.7) 

6f,i = -(?-3) (18.8) 

b F ,2 = q-3 (18.9) 

and 

&F )3 = -(<Z-2)(g-3) . (18.10) 
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FIG. 15: Graph F 7 . 



The chromatic polynomial may also be expressed in the form of Eq. (2.9), with j max = 3, 
namely 

3 

P(F m ,q) = y £ i c F j(\ F j) m . (18.11) 

Using Eq (2.13) (or (2.15)) of Ref. one can calculate the coefficients c F j for j = 1,2,3 
from the generating function. Specifically, we have 

(a FQ X F1 + a F iX F1 + a F2 ) , 1Q1 „, 

CF ' X = 7a — r~v~\ — 3rr (l8 - 12) 

{A F i — A F 2)\A F i — A F 3) 

(a F0 X 2 F2 + a F iX F2 + a F2 ) , 1D1 „, 

CF - 2 = 7a — r~v~\ — 3rr (l8 - l3) 

(,AF,2 — A F ,l){A F ,2 — AF,3) 

and 

{a F0 X F3 + a F1 X F3 + a F2 ) 

CF ' 3 = 7a — r~v\ — 3rr • (l8 - l4) 

\A F ,3 — A F i)[A F 3 — Af,2) 

As discussed before [HHTU], the X F jS appear via the factorized form of the denominator of 
the generating function, 

3 

V(F,q,x) = Y[(l - X Fj x) . (18. 15) 

3=1 

Equivalently, the X F jS are determined from the equation £ 3 + b F ^ 2 + b F2 C, + b F3 = 0, i.e., 

e + (3 - q)e + (q - 3)C - (q - 2){q - 3) = . (18.16) 
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Let us define 



and 



R F = 3(4g 3 - 24g 2 + 76g - 93) 



1/3 



'18.17) 



'18.18) 



S F = \A(q - 3) \2q 2 + 6q - 9 + 3y/R F 
With appropriate choices of branch cuts for the various fractional powers in (18.18), we have 

2(q 



A 



S f 

~6~ 



-3)(g 
35* f 



18.19) 



The other Ai?j, j = 2, 3 can be written explicitly in a similar manner. Thus, this family 
is valuable as an illustration of a family of planar triangulation graphs with a chromatic 



polynomial of the form (2.9) and with A terms that are different from those in (3.1) and 



indeed, are nonpolynomial, in contrast to the families with chromatic polynomials of the 



form (2.8) or (3.16). 



With regard to the evaluation of P(F m , q) at q = r + 1 (with an appropriate choice of 



branch cuts for the square and cube roots), A^i and one of the other two roots of (18.16) 
comprise the complex-conjugate pair 



- 1 + V5 ± (-38 + 18X/5) 1 



/2, 



;i8.20) 



with magnitude 0.485867.., while the third root of (18.16) is equal to —1. Since 0.485867.. 
is less than r — 1 = 0.6180.., the corresponding two coefficients do not have to, and do not, 
vanish at q = r + 1. Since the third root has magnitude greater than r — 1, its coefficient 
must vanish at q = r + 1 in order for \P(F m ,r + 1)| to obey the Tutte upper bound (1.1). 
With these values of the Xf,/s at q — r + 1, the ratio r(F m ) vanishes (exponentially rapidly) 
as m — > oo and r(F DO ) = 0. This illustrates the general property that if G ptjm is a family of 
planar triangulations with P(G pttm ,q) of the form (2.9) and a = 1 in (2.1), and if none of 



the Xa pt ,j has magnitude equal to r — 1 when evaluated at q 



r + 1, then, since (i) the Ac 



with |Ag t j\ > t— 1 have coefficients that must vanish, and (ii) the Xg^j with \Xq, 



P t ,3 



' P t,3 



P t,3 

< r-1 



give zero contribution in the limit m — > oo, it follows that r(G pt ) = 0. We calculate 

op = 0.786151.. 



18.21' 



The term A^i is real and positive and is dominant for q > r + 2 = 3.618.... In this 
interval, the other two roots, Xpj, j = 2,3 are complex, with smaller magnitudes. At 
q = t + 2, Xp,i = —1 and |Ai? 2 | = | A^ ;3 | = 1, so all Xp,j are degenerate in magnitude. Hence, 
in the notation of [S], q c — r + 2 for this family. At q = 3, all Xpj = 0, j = 1,2,3, as is 



obvious from Eq. (18.16). 
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We exhibit the first few P(F m ,q). For m = 1, P(Fi,q) = clf,o, as given above in (18.5). 
For m = 2 to m = 6, 

P(F 2 , q) = q(q -l)(q- 2)(q 3 - 9q 2 + 29q - 32) (18.22) 

P{F 3 , q) = q{q - l){q - 2){q - 3)(g 3 - 9q 2 + 30q - 35) (18.23) 
P(F 4 , q) = q(q - l)(q - 2)(q - 3)(g 4 - 12g 3 + 58g 2 - 133g + 119) (18.24) 

P(F 5 , q) = q(q - l)(q - 2)(q - 3)(g 5 - 15g 4 + 95g 3 - 317g 2 + 553g - 398) 

(18.25) 

P(F 6 , q) = q(q - l)(q - 2)(q - 3) 2 (g 5 - 15g 4 + 96g 3 - 327g 2 + 591g - 447) 

(18.26) 

As m increases further, P(F m , q) has increasingly high powers of the factor (q — 3). 

As with the other planar triangulation families, the F m family has chromatic zeros near 
to t + 1. We find that these approach r + 1 as m gets large. Depending on the value of m, 
P(F m , q) also may have real zeros in the interval [q w , 3). The complex zeros of P(F m , q) form 
a complex-conjugate arc, with arc endpoints at the complex zeros of Ftp, namely q, q* ~ 
1.9111 ±2.65021 



XIX. SOME IMPLICATIONS FOR STATISTICAL PHYSICS 



One of the interesting aspects of the present work is its implications for nonzero ground 
state entropy of the Potts antiferro magnet. (For background on the Potts model, see Refs. 
[35| 136] . [8].) This stems from the identity noted above, P(G,q) = Z(G,q,T = 0)paf = 
Wtot(G,q). As above, we denote the formal limit of a family of graphs G as n(G) — > oo by 
the symbol {G}. We recall that the entropy per vertex is given by So = kp In W, where W 
is the degeneracy per vertex, related to the total degeneracy of spin configurations of the 
zero-temperature Potts antiferromagnet (or equivalently the number of proper g-colorings 
of the graph) by W tot by W({G},q) = lim n _ >00 [H 7 " tot (G, q)] l l n . We refer the reader to Ref. 
[S] for a discussion of a subtlety in this definition resulting from a certain noncommuativity 
that occurs for a special set of values of q, denoted as {q s }, namely 

lim lim [P{G,q)] 1/n ^ lim lim [P(G, q)] 1/n . (19.1) 

q— >q s n— >oo n— >oo q— >q B 
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For the one-parameter families of planar triangulations considered here, this set of special 
values {q s } includes q = 0, 1, 2, r + 1 and, for cases where the chromatic number is 4, also 
q = 3. Because of this noncommutativity, it is necessary to specify the order of limits taken 
in defining W. For a particular value q = q s , we thus define 



]l/n 

" qny J ) HS) — [ J K^iH) 

q— >q a n— yoo 

and 



W qn ({G}, q s ) = lim lim [P{G, q)] l ' n (19.2) 

9— >9s fi— >oo 



W nq ({G}, q s ) = lim lim [P(G, q)] 1/n . (19.3) 

For real g > x(Gpt )m ), both of these definitions are equivalent, and in this case we shall write 
W qn ({G},q s ) = W ng ({G},q s ) EE W({G},q s ). 

We generalize our calculations in [1] as follows. First, for a family of planar triangulations 



Gpt, m with P(G p t :rn ,q) having the form (2.9) with j max — 1, W — [AgjJ 1 ■ As an example, 
consider the family 

R m = P m + P 2 , (19.4) 

with n = m + 2 considered in [I]. Here P m is the path graph with m vertices and the join 
G + H of two graphs G = (Vg, Eg) and H = (Vh, Eh) is defined as the graph with vertex 
set Vg+h = Vq U Vh and edge set Eq+h comprised of the union of Eq U -Eh with the set 
of edges obtained by connecting each vertex of G with each vertex of H . Thus, R\ = K3, 
i?2 = K4, etc. The W function for the m — > 00 limit of this family is W({R},q) = q — 3, 
and S > for q > 4. For the family TC m 

W({TC},q) = (\ TC ) 1/3 . (19.5) 



where \tc was given in (3.24). The function Arc is a monotonically increasing function of q, 
which passes through zero at q — 2.54660... and increases through unity at q = 3 so that (for 
the m — > 00 limit of this family) Sq > for q > 3. For the family I m of iterated icosahedra, 



where 



W({I},q) = (X I ) 1 /\ (19.6) 

\j = ( q - 3)(g 8 - 24g 7 + 260g 6 - 1670g 5 + 6999g 4 - 19698g 3 

+ 36408g 2 - 40240g + 20170) . (19.7) 

The function Xj vanishes at three real values of q, namely q = 2.618197.. (i.e., slightly above 
t + 1), at q — 3, and at q = 3.222458... This function is positive for q > 3.222458.. (as well 
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as in an interval 2.618.. < q < 3) and increases through unity as q increases through the 
value 3.5133658.., so that in this latter interval, So > 0. 

A second general result is that for a family of planar triangulations G pt ^ m with P(G pt>m , q) 
having the form (3.16), it follows that (i) W qn ({G pt },q) = q — 2 for q > 3; (ii) even in the 
presence of noncommutativty, W({G pt }, q) = q — 2 for q > 4, so that S > in this interval 
(and also in the interval q > 3 if one uses W qn ({G pt }, q)). This result applies, in particular, to 
the families B m , H m , L m , -D m _4,2, and _D m _4 5 3. Although the family P(D m ^ m , q) is of the form 
( 2.9[ ) with j max = 6, the dominant term for q > x{D m>m ) is again q — 2, so that in this interval 
W({Dd},q) = q — 2 for this family also. In contrast, P(S mit m 2 ,q) has cs,u = and hence 
lacks the term that would normally be dominant as m,\ or rri2 goes to infinity. In this case, for 
q > 4 where there is no noncommutativity in limits, we find W({S},q) = yV(q — 2)(q — 3) . 

For the family F m , we find that Ai in Eq. (18.19) is dominant for q > q c = r + 2 = 3.618..., 
so that in this region, 

W({F},q) = \ Ffl (19.8) 

Furthermore, since Af,i > 1 for real q > r + 2, it follows that So > for (the m — > oo limit 
of this family of graphs) for this range q > r + 2. 

For a regular lattice graph G it is of interest to investigate the dependence of W({G}, q) 
on the vertex degree (coordination number) d. This study was carried out in [S] 13"8"] . 
and it was shown that W({G},q) is a non-increasing function of d. This is understood as 
being a consequence of the fact that (except for tree graphs, which are not relevant here), 
roughly speaking, increasing the vertex degree tends to increase the constraints on a proper 
g-coloring of the vertices and therefore tends to decrease W({G},q). One is also motivated 
to investigate the same question with the families of planar triangulations under study here. 
However, since d e ff = 6 for a family of planar triangulation graphs one is limited to a fixed 
d e ff = 6 and hence cannot carry out the type of comparative study involving a variation in 
d e ff that was performed in [S], [3H EH]- In [I] and the present work, we have found that 
families of planar triangulations can have different IV({G}, q) functions. This is consistent 
with the results in [8j-[38j. Indeed, one has already encountered examples of this. For 
example, the square and kagome (3 ■ 6 • 3 ■ 6) lattices both have the same vertex degree, 
namely 4, but they have different W functions, and similarly, the honeycomb, (3 • 12 2 ), and 
(4 • 8 2 ) lattices have the same vertex degree, namely 3, but again, they have different W 
functions PE3EH]. 
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XX. COMPARATIVE DISCUSSION 



In this section we give a comparative discussion of some limiting quantities for the various 
families of planar triangulations that we have studied so far. For one-parameter families of 
planar triangulations G pt ^ m for which P(G pt>m ,q) is of the form (2.8) we have proved that 
r {G p t,oo) — and have investigated the various values of ac pt defined in (3.18). This constant 
is strictly less than unity, and it is of interest to see which families yield larger and smaller 
values of ac vt - We display the values that we have obtained in Table VII As is evident, in 
the set of j max = 1 families of planar triangulations, the family of cylindrical strips of the 
triangular lattice, TC m (equivalently, iterated octahedra) yields the largest value of ac pt , 
which is within 9 % of its upper bound of 1. In the j max = 3 families, the one that yields 
the largest value of the limiting ratio r(G p4j00 ) is the family, B m , with r^B^) = 0.6180.. 
A second type of asymptotic limiting function is W({G},q). We have given a comparative 
analysis of this in the previous section. 

We have also investigated the values of P(G p t, m , q) at q = x(G p t,m) for the families of 
planar triangulations that we have studied. Recall the definition that a graph G is fc-critical 
iff x{G) = k and P(G, k) = k\. We find a variety of behavior. For example, (i) x(i? m ) = 4 
and P(i? m ,4) = 4!, so R m is 4-critical; (ii) x(TC m ) = 3 and P(TC m ,3) = 3!, so TC m is 

3- critical; but (iii) x = 4 for I m , H m , L m , -D m _ 4 3 , and F m > 3 but none of these families is 

4- critical. For other families G pt ,mi the chromatic number depends on whether m is even or 
odd. For example, for even m, \ — 3 for B m and -D m _ 4 2 , and these graphs are 3-critical, 
while for odd m, x = 4 for B m and -D m _ 4i2 , but neither of these graphs is 4-critical. 

There are a number of further directions in which the present research could be extended. 
One could, for example, study Tutte polynomials, or the equivalent, Potts model partition 
functions Z(G pt ,q,v), for planar triangulation graphs G pt (where v is a temperature-like 
Boltzmann variable). However, a number of the special properties that make the chromatic 
polynomials of these graphs amenable to analysis do not generalize to the full Potts model 
partition function. These include, for example, the use of the complete-graph intersection 
theorem and the property that the chromatiic polynomial has q(q — l)(q — 2) as a factor. For 
example, for the triangle graph K3 itself, the Potts model partition function Z(K 3 , q, v) = 
(q + v) 3 + (q — l)v 3 = q(q 2 + 3qv + 3v 2 + v 3 ) only has the factor q. Another direction of 
investigation would be to calculate weighted-set chromatic polynomials [39l HO] and Potts 
model partition functions with an external magnetic field that favors or disfavors a single 
value of q or a set of s such values [4T] for planar triangulation graphs. Work on this is 
underway. 
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TABLE VII: Some asymptotic limiting quantities for one-parameter families of planar triangulations. The 
shorthand notation 3me, 4mo means \ = 3 if m is even and x = 4 if m is odd. Additional information for x 
values is x(^3) = 3, x(-^i) = 4, and x(-^2) = 3. Numerical values are quoted to three significant figures. 





n '\S'pt,m) 




Jmax 




a G pt 




m + 2 


4 if m > 2 


1 





(—1 + v5)/2 = 0.618 




3m 


3 


1 





(3 — V5 ) ' = 0.914 


Ira 


9m + 3 


4 


1 





[(—315 + 141V5 )/2J 7 = 0.8055 


GcM,m 


8m + 3 


3 


1 





[(115 - 51V5 )/2] 1/s = 0.885 




i /I 

m + 4 


A " J? \ O 

4 if m > 3 


3 





0.786 




m + 2 


3me, 4mo 


3 


(-1 + V5 )/2 = 0.618 


1 


IT 


m + 5 


4 


Q 

o 


(/ — ovO )l £ — U.14D 


i 
i 


L m 


m + 5 


4 


3 


-2 + ^5 = 0.236 


1 




m + 5 


3me, 4mo 


3 


-4 + 2^/5 = 0.472 


1 


An-4,1 


m + 6 


4 


3 


(-15 + 7^5 )/2 = 0.326 


1 


A),m-2 


m + 7 


3me, 4mo 


3 


-13 + 6\/5 = 0.416 


1 


-t>l,m-2 


m + 8 


4 


3 


-22 + 10V5 = 0.361 


1 


D m ,m 


2m + 9 


3me, 4mo 


6 


(3-V5 )/2 = 0.382 


1 


Sm,m 


m + 7 


4 


4 





(-1 + \/5 )/2 = 0.618 



XXI. CONCLUSIONS 

In this paper we have presented an analysis of the structure and properties of chro- 
matic polynomials P(G pt ^, q) of families of planar triangulation graphs G ptt rn, where 
m = (m 1 ,...,m p ) is a vector of integer parameters. We have discussed a number of spe- 
cific families with p = 1 and p = 2. These planar triangulation graphs form a particularly 
attractive class of graphs for the analysis of chromatic polynomials because of their special 
properties. One of these is the fact that when evaluated at q — r + 1, the chromatic polyno- 



mial of a planar triangulation graph satisfies the Tutte upper bound (1.1). We have studied 



the ratio of |P(G p t )rf i, r + 1) | to the Tutte upper bound (r — 1)™ 5 and have calculated limiting 
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values of this ratio as n — > oo for various families of planar triangulations. We also have used 
our calculations to study zeros of these chromatic polynomials. Among our results, we have 
shown that if G vt ^ is a planar triangulation graph with a chromatic polynomial P(G pt ,rn, q) 
of the form (8.1 ), then (i) the coefficients c q % must satisfy a number of properties, which we 
have derived; and (ii) P(G p t^, q) has a real chromatic zero that approaches (1/2) (3 + y/5 ) 
as one or more rrii — > oo. We have constructed a p = 1 family of planar triangulations 
with real zeros that approach 3 from below as m — > oo. A one-parameter family F m with 
jmax = 3 and nonpolynomial Xp,j has been studied. We have also presented results for 
a number of results for chromatic polynomials of various two-parameter families of planar 
triangulations. Implications for the ground-state entropy of the Potts antiferromagnet are 
discussed. Our results are of interest both from the point of view of mathematical graph 
theory and statistical physics and further show the fruitful connections between these fields. 
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